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t  distribution  of  radiation  from  that  surface.  The  existence  of  such  a  transfer 
function  implies  a  shift-invariant  scattering  function  which  does  not  change 
shape  with  the  angle  of  the  incident  beam.  This  result  greatly  reduces  the 
quantity  of  data  required  to  completely  characterize  the  scattering  properties 
of  a  surface.  For  a  large  class  of  well-behaved  surfaces  this  transfer  func¬ 
tion  is  described  in  terms  of  only  the  rms  surface  roughness  and  the  surface 
autocovariance  function.  It  thus  provides  a  straightforward  solution  to  the 
inverse  scattering  problem  (i.e.,  determining  surface  characteristics  from 
scattered  light  measurements) .  Once  the  surface  characteristics  are  known,  the 
same  theory  provides  an  equally  simple  method  of  predicting  the  wavelength 
dependence  of  the  scattered  light  distribution. 

An  extensive  experimental  program  has  accompanied  this  theoretical  develop¬ 
ment.  <  The  apparatus  and  experimental  procedures  utilized  in  measuring  the 
angular  distribution  of  light  scattered  from  a  variety  of  optical  surfaces 
for  several  different  angles  of  incidence. and  wavelengths  are  described  in 
detail.  Experimental  verification  of  the  shift-invariant  scattering  function 
has  been  successfully  demonstrated  for  smooth  surfaces  (oj/«X).  The  scattered 
light  measurements  from  rough  (diffusely  reflecting)  surfaces  results  in  a 
scattering  function  which  is  shift-invariant  over  only  a  small  range  of  angles 
and  departs  significantly  from  the  predicted  behavior  at  large  scattering 
angles. 

A  computer  program  has  been  developed  that  operates  upon  scattered  light 
data  to  yield  the  total  integrated  scatter,  the  surface  transfer  function,  the 
rms  surface  roughness,  and  the  surface  autocovariance  function.  Although 
accurate  determination  of  microstructure  on  optical  surfaces  is  extremely  dif¬ 
ficult  to  accomplish  by  direct  measurement  (thus  the  motivation  for  attempting 
to  solve  the  inverse  scattering  problem),  favorable  comparisons  of  predicted 
surface  characteristics  with  the  corresponding  measured  quantities  have  been 
demonstrated  for  both  smooth  surfaces  and  moderately  rough  surfaces.  In 
addition,  experimental  verification  of  the  inverse  scattering  program  was 
accomplished  indirectly  by  supplying  scattered  light  data  of  one  wavelength  ns 
input  to  the  inverse  scattering  program  in  order  to  determine  the  relevant 
surface  characteristics;  then  this  information  was  used  to  predict  the  scat¬ 
tering  function  at  a  different  wavelength.  Excellent  agreement  with  the 
measured  scattering  function  at  that  wavelength  was  achieved. 

Since  the  above  technique  involves  numerical  computations  on  sampled  data, 
an  analytical  expression  for  a  wavelength  scaling  law  is  not  required  to  deter¬ 
mine  the  scattering  function  at  any  desired  wavelength.  However,  in  order  to 
gain  insight  into  the  wavelength  dependence  of  surface  scatter  phenomena,  a 
wavelength  scaling  law  for  smooth  surfaces  was  derived  and  verified.  This 
scaling  law  consists  of  a  change  in  the  scattering  angle  as  well  as  a  change 
in  the  amplitude  of  the  scattering  function  with  changes  in  wavelength.  It 
therefore  provides  a  valuable  tool  for  predicting  the  scattering  behavior  in 
certain  angular  regions  or  wavelength  ranges  where  direct  measurements  are 
difficult  to  obtain. 
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CHAPTER  1 

INTRODUCTION 

The  relationship  between  surface  micro-structure  and  radiant 
energy  scattering  plays  an  important  role  in  many  areas  of  technical 
interest.  These  include  the  trade-off  between  cost  and  performance  in 
the  fabrication  of  optical  surfaces,  design  considerations  for  stray- 
light  rejection  systems,  evaluation  of  machined  metal  mirrors  for  high- 
energy  laser  applications,  laser-radar  backscatter  signature  programs, 
and  a  host  of  other  applications  requiring  extensive  scattering  data. 

If  the  scattering  mechanism  were  completely  understood,  surface  prepa¬ 
ration  techniques  or  measurement  programs  in  many  of  these  areas  could 
possibly  be  changed  to  obtain  more  favorable  results. 

Background 

If  a  propagating  wave  is  incident  upon  a  perfectly  plane  sur¬ 
face,  the  reflected  wave  is  concentrated  in  the  specular  direction  as 
determined  by  the  well-known  laws  of  reflection.  Another  idealized 
surface  is  the  perfectly  diffuse  reflector  which  scatters  light  accord¬ 
ing  to  Lambert’s  cosine  law.  A  more  physically  realistic  situation  is 
shown  in  Fig.  1.1,  which  illustrates  the  optical  scattering  that  occurs 
when  light* is  reflected  from  a  rough  surface.  If  the  surface  is  not 
too  roughv'the  reflected  light  consists  of  a  specular  component  plus  a 
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SPECULAR  BEAM 


Fig.  1.1.  Schematic  Representation  of  Reflectance 
from  a  Rough  Surface. 


diffuse  component  which  is  scattered  over  a  wide  range  of  angles  cen¬ 
tered  upon  the  specular  beam. 

One  of  the  earliest  investigators  of  scattering  from  a  rough 
surface  was  Lord  Rayleigh.  In  1896  (Rayleigh,  1945)  he  was  investiga¬ 
ting  the  reflection  of  acoustic  waves,  and  later  (Rayleigh,  1901)  he 
noted  the  effects  of  poorly  polished  surfaces  on  optical  performance. 

He  examined  the  effects  of  surface  roughness,  wavelength,  and  angle  of 
incidence  on  the  reflected  beam,  Chenmoganadam  (1919)  derived  a  theory 
of  scattered  light  based  on  the  phase  shift  of  the  reflected  beam  due 
to  the  rough  surface. 
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However,  it  was  not  until  the  problem  of  background  clutter  in 
radar  applications  became  apparent  that  a  determined  effort  was  made 
to  solve  the  scattering  problem  for  random  surfaces.  For  example,  scat¬ 
tering  from  the  sea  motivated  the  work  of  Davies  (1954)  as  well  as 
others  (ilakc,  1950;  Barrick,  1970;  Bass,  1968;  Beard,  1961;  and  Fuks, 
1966).  Considerable  work  has  also  been  done  in  attempts  to  explain 
radar  reflection  from  the  moon  (Daniels,  1961;  Evans  and  P&ttengill, 
1963;  Fung  and  Moore,  1964;  Fung,  1967;  and  Hagfors,  1964). 

Random  rough  surfaces  have  been  treated  in  two  different  ways. 
Rough  surfaces  made  up  of  a  random  array  of  objects  or  shapes  with 
known  scattering  characteristics  were  investigated  by  Ament  (1960) , 
Twersky  (1957),  Spetner  (1958),  and  Peake  (1959).  The  other  approach 
taken  by  Isakovich  (1952),  Ament  (1953),  Eckart  (1953),  Feinstein 
(1954),  fiavies  (1954),  and  Beckmann  (1957),  treats  the  rough  surface  as 
a  stochastic  process. 

S.nce  optical  surfaces  clearly  fall  into  the  second  classifica¬ 
tion  of  random  surfaces,  Bennett  and  Porteus  (1961)  expanded  and  experi¬ 
mentally  investigated  the  scattering  theory  of  Davies  (1954) .  From  this 
and  subsequent  work  (Bennett,  1963;  and  Porteus,  1963)  the  reflectance 
properties  of  samples  with  a  measured  surface  roughness  were  directly 
compared  to  theory  with  good  results.  Interest  in  these  measurements 
led  to  investigations  at  the  Optical  Sciences  Center,  University  of 
Arizona,  by  Mott  (1971),  McKenney,  Orme  and  Mott  (1972),  Orme  (1972), 
DeBell  and  Harvey  (1974),  Shack  and  DeBell  (1974),  and  Shack  and 
Harvey  (1975). 
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The  Bidirectional  Reflectance  Distribution  Function  (BRDF)  was 
introduced  by  Nicodemus  (1970)  as  a  quantity  which  completely  describes 
the  reflectance  (or  scattering)  properties  of  a  given  surface. 
Bidirectional  reflectance  data  for  spectral  regions  extending  into  the 
infrared  have  been  collected  in  connection  with  heat-transfer  analysis. 
Another  area  of  interest  involves  the  possibilities  for  spectro- 
chemical  analysis  by  reflected  radiation,  principally  in  connection 
with  remote  sensing  of  the  earth  and  other  planetary  surfaces.  Also 
BRDF  measurements  have  been  made  on  many  proposed  baffle  materials  for 
use  in  the  Large  Space  Telescope  program  (Breault  and  Fannin,  1973). 
Recently  there  has  been  a  great  deal  of  activity  in  the  area  of  BRDF 
measurements  of  machined  metal  mirrors  to  be  used  in  high-energy  laser 
applications  (Young,  1975;  Curcio,  1975;  Decker,  Bennett,  and  Bennett, 
1975;  Church  and  Zavada,  1975;  and  Stover,  1975). 

Dissertation  Content 

In  this  dissertation  the  scattering  of  light  from  optical 
elements  is  considered  to  be  solely  a  surface  phenomenon.  Light¬ 
scattering  from  optical  surfaces  is  then  treated  as  a  diffraction 
process  in  which  the  pupil  function  has  random  phase  variations  in 
addition  to  any  existing  amplitude  variations.  A  complete  Fourier 
treatment  of  near-field  scalar  diffraction  theory  is  therefore  devel¬ 
oped  in  Chapter  2.  This  diffraction  theory  is  generalized  in  Chapter  3 
to  include  phase  perturbations  that  lead  to  scattered  radiation.  Appro 
priate  assumptions  are  then  made  concerning  the  statistical  properties 
of  optical  surfaces  and  an  analytical  expression  is  obtained  for  the 
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transfer  function  of  a  scattering  surface.  The  existence  of  such  a 
transfer  function  implies  a  shift-invariant  scattering  function  which 
does  not  change  shape  with  the  angle  of  the  incident  beam.  This  result 
greatly  reduces  the  quantity  of  data  required  to  completely  character¬ 
ize  the  scattering  properties  of  a  surface.  For  a  large  class  of  well- 
behaved  surfaces  this  transfer  function  is  described  in  terms  ot  only 
the  rms  surface  roughness  and  the  surface  autocovai ranee  function.  It 
thus  provides  a  straightforward  solution  to  the  inverse  scattering 
problem  (i.e.,  determining  surface  characteristics  from  scattered  light 
measurements).  Once  the  surface  characteristics  are  known,  the  same 
theory  provides  an  equally  simple  method  of  predicting  the  wavelength 
dependence  of  the  scattered  light  distribution. 

An  extensive  experimental  program  has  accompanied  this  theoreti¬ 
cal  development.  Chapter  4  describes  in  detail  the  apparatus  and  experi¬ 
mental  procedures  utilized  in  measuring  the  angular  distribution  of  lignt 
scattered  from  a  variety  of  optical  surfaces  for  several  different  angles 
of  incidence  and  wavelengths.  The  results  of  these  experiments  are 
reported  in  Chapter  S  and  compared  to  theoretical  predictions  with 
generally  good  agreement. 

Consistent  with  most  research  efforts,  not  all  questions  con¬ 
cerning  this  topic  are  completely  answered  in  this  dissertation.  After 
a  brief  summary  of  results,  Chapter  6  is  therefore  devoted  to  a  few 
comments  concerning  new  theoretical  considerations  and  suggestions  for 
future  research. 


CHAPTER  2 


A  FOURIER  TREAWENT  OF  NEAR-FIELD 
SCALAR  DIFFRACTION  THEORY 

N 

The  phenomenon  of  diffraction  involves  a  wave  field  incident 
upon  one  or  more  objects  or  apertures  with  absorbing  or  conducting 
surfaces.  The  calculation  of  the  wave  field  emerging  from  such  a 
diffracting  system  is  the  goal  of  all  diffraction  theories. 

It  should  be  emphasized  that  both  the  Kirchhoff  and  Rayleigh- 
Sommerfeld  theories,  as  well  as  the  present  discussion  in  this  paper, 
treat  light  as  a  scalar  phenomenon.  (For  a  detailed  treatment  cf  the 
historical  development  of  diffraction  theory,  see  Goodman,  1968, 
pp.  30-S6.)  Such  an  approach  entirely  neglects  the  fact  that  the 
various  components  of  the  electric  and  magnetic  field  vectors  are 
;oupiea  through  Maxwell’s  equations  and  cannot  be  treated  independently. 

Microwave  experiments  have  shown  that  scalar  theory  yields 
very  accurate  results  provided  that:  (1)  the  diffracting  aperture  is 
large  compared  to  a  wavelength,  and  (2)  the  diffracted  wave  field  is 
observed  far  from  the  aperture.  It  is  significant  that  although  the 
present  treatment  is  limited  by  being  a  scalar  theory,  the  above 
approximations  are  not  imposed  during  the  mathematical  formulation  as 
they  are  in  the  Kirchhoff  theory.  Furthermore,  the  following  devel¬ 
opment  provides  much  more  insight  than  the  conventional  Rayleigh- 
Sommerfeld  theory. 


The  Diffracted  Wave  Field 
as  a  Superposition  of  Plane  Waves 


The  fundamental  diffraction  problem  consists  of  two  parts:  (1) 
determining  the  effect  of  introducing  the  diffracting  screen  upon  the 
field  immediately  behind  the  diffracting  screen,  and  (2)  determining 
how  it  affects  the  field  downstream  from  the  diffracting  screen  (i.e., 
what  is  the  field  immediately  behind  the  diffracting  screen  and  how 
does  it  propagate) . 

Consider  first  the  propagation  problem  and  let  the  complex 
amplitude  distribution  of  the  optical  disturbance  in  plane  PQ  be  repre¬ 
sented  by  the  scalar  function  0).  This  scalar  disturbance  in  PQ 

will  be  considered  the  only  radiation  contributing  to  the  field 
£/(£,$;£)  in  plane  P  (see  Fig.  2.1),  2  has  a  parametric  relationship 
since  it  is  a  function  of  the  observation  plane.  Note  that  a  scaled 
coordinate  system  is  utilized  in  which  &  -  x/X,  p  =  y/X,  z  =  z/X. 

Initial  Conditions 

It  will  be  assumed  that  the  complex  amplitude  of  any  monochro¬ 
matic  optical  disturbance  propagating  through  free  space  must  obey  the 
time- independent  wave  equation  (Helmholtz  equation).  We  will  also 
assume  that  the  Fourier  transform  of  the  scalar  field  £/o(A,£;0)  exists. 
This  is  not  a  severe  restriction,  however,  as  Bracewell  (1965)  points 


out  that  physical  possibility  is  a  valid  sufficient  condition  for  the 
existence  of  a  Fourier  transform. 


The  Direct  Application  of  Fourier  Transform  Theory 

We  can  thus  define  the  following  Fourier  transform  relation- 

» 

ships  that  exist  for  planes  Pq  and  P. 

00 

V«.e;0)  B  ||  vo(&,9i0)  e'^^^d&dg  (i) 

-00 

00 

0)  s  ||  et'2n(aA+e^)dadB  (2) 

-00 

00 

A(a,S;&)  *  ||  Vi&Jih)  e-x2Ha*+Wd&d9  (3) 

-00 

00 

V(x,ir,z)  =  ||  A( a,M)  ei2v dad& .  (4) 

—  00 

Equations  (2)  and  (4)  indicate  that  the  monochromatic  scalar 
wave  field  in  planes  Pq  and  P  can  be  decomposed  into  plane  wave  compo¬ 
nents  whose  amplitudes  are  a  function  of  the  direction  cosines  of  the 
propagation  vector.  The  functions  Ao(o,6;0)  and  A (a, 6; 2)  will  be 
referred  to  as  the  direction  cosine  spectrum  of  plane  waves  contributing 
to  the  disturbance  U0(&,9', 0)  and  U(£,p;2)  respectively.  The  direction 


cosine  spectrum  of  plane  waves  is  used  here  in  lieu  of  the  angular 
spectrum  of  plane  waves  discussed  by  Ratcliffe  (1956)  and  others.  This 
is  consistent  with  a  more  general  treatment  which  is  not  restricted  to 
small  angles. 

In  the  scaled  coordinate  system  V2  «  X2V2,  and  k2  »  X2fc2  • 

(2tt) 2 .  Hence  the  Helmholtz  equation  becomes 

[V2  ♦  (2ir)2]tf (A,#;*)  «  0. 


(5) 


Now  by  applying  Eq.  (4)  and  requiring  the  individual  plane  wave 
components  to  satisfy  the  Helmholtz  equation,  we  find 


where 


4(a,0;&)  ■  4o(a,0;O)  e 


/l 


i2vy& 


(6) 


The  Transfei  Function  of  Free  Space 

Since  Eq.  (6)  relates  the  Fourier  transforms  of  the  scalar 
fields  in  planes  Pq  and  P  it  can  be  rewritten  in  terms  of  a  transfer 
function  for  free  space,  tf(a,0;S) 

s  *  e<2’YS-  <7> 
We  have  thus  far  applied  no  restrictions  on  y  and  two  regions 


of  interest  are  apparent:  that  for  real  values  of  y  and  that  for 
imaginary  values, 


y  =  n  -  (a2+g^) 


for  (a2  +  32)  si  y  is  real 


for  (a2  +  32)  >  1 


Consider  now  a  unit  circle  in  the  a- 3  plane  of  direction  cosine 
space  as  shown  in  Fig.  2.2.  Inside  this  unit  circle  y  is  real  and  the 
corresponding  part  of  the  disturbance  will  propagate  and  contribute  to 
the  wave  field  in  plane  P.  However,  those  components  of  the  direction 
cosine  spectrum  which  lie  outside  the  unit  circle  have  imaginary  values 
of  y  and  represent  that  part  of  the  disturbance  which  experiences  a 
rapid  exponential  decay.  This  is  the  part  of  the  disturbance  which  is 
commonly  referred  to  as  the  evanescent  wave  (Goodman,  1968) . 

Let  U0(x,y; 0)  be  the  product  of  the  complex  amplitude  transmit¬ 
tance  of  a  diffracting  screen  and  the  complex  amplitude  distribution 
incident  upon  it.  Figure  2.3(a)  illustrates  this  quantity  broken  down 
into  the  part  which  propagates  and  the  part  which  makes  up  the  evanes¬ 
cent  wave  for  the  case  of  a  unit  amplitude  plane  wave  normally  incident 
upon  a  circular  aperture.  The  direction  cosine  spectrum  of  plane  waves 
associated  with  these  respective  optical  disturbances  are  shown  in  Fig. 


2.3(b). 


Note  that  the  sharp  corners  on  the  original  disturbance  in 


Fig.  2.3(a)  correspond  to  Kirchhoff's  unnecessary  boundary  conditions. 
It  is  the  propagating  part  only  that  accurately  represents  the 
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The  plane  wave  components  msi 
circle  will  propagate,  and  the 
wave  components  outside  this  c 
contribute  to  the  evanescent  w 
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disturbance  immediately  behind  the  diffracting  aperture  which  will  con¬ 
tribute  to  the  disturbance  downstream. 

It  is  now  clear  that  the  complex  amplitude  distribution  in  plane 
P  can  be  determined  by  Fourier  transforming  the  original  disturbance 
yo(£»p;0),  then  multiplying  the  resulting  direction  cosine  spectrum  of 
plane  waves  4o(a,8;0)  by  the  transfer  function  of  free  space  given  in 
Eq.  (7),  and  finally  by  applying  the  inverse  Fourier  transform  integral 
of  Eq.  (4).  However,  the  limits  of  integration  on  Eq.  (4)  must  be 
changed  such  that  the  integration  is  performed  only  over  the  unit 
circle  instead  of  over  the  entire  a-B  plane. 

The  above  analysis,  in  which  an  optical  disturbance  is  repre¬ 
sented  as  a  superposition  of  plane  waves,  corresponds  to  the  transfer 
function  approach  in  image  formation  and  yields  considerable  insight 
into  the  behavior  of  these  plane  wave  components  during  the  phenomenon 
of  diffraction. 


The  Diffracted  Wave  Field 
as  a  Superposition  of  Spherical  Waves 

The  convolution  theorem  (Bracewell,  1965)  of  Fourier  transform 

theory  requires  that  a  convolution  operation  exists  in  the  domain  of 


real  space  that  is  equivalent  to  Eq.  (6) 


The  Point  Spread  Function 

We  thus  have  the  alternative  method  of  expressing  the  complex 
amplitude  distribution  in  the  observation  plane  by  the  convolution  of 


# 

the  original  disturbance  with  a  point  spread  function.  The  point 
spread  function  is  obtained  by  taking  the  inverse  Fourier  transform  of 
the  transfer  function  found  in  Eq.  (7). 

Starting  with  the  well-known  Weyl  expansion  formula  (Weyl,  1919), 
Lalor  (1968)  obtained  a  result  which,  with  straightforward  modification, 
yields 


j  ei2lTY*  e<2,(a4+^<kU 


_1_ 

2tt 


i2  Ttf' 
e 


(9) 


where 

£2  =  A2  +  p2  +  S2. 

The  left  side  of  Eq.  (9)  is  the  inverse  Fourier  transform  of  the  trans¬ 
fer  function  of  free  space.  The  appropriate  point  spread  function  is 
thus  given  by 

«*.»!*)  ■  *■-'  {ei2'Y8}  ■  (s?  -  <)  £  *1^  •  »°> 


Huygens'  Principle 

Recall  now  the  assertion  by  Christiaan  Huygens  (Thompson,  1912) 
in  1678  that  each  element  of  a  wavefront  may  be  regarded  as  the  center 
of  a  secondary  disturbance  which  gives  rise  to  spherical  wavelets;  and 
moreover  that  the  position  of  the  wavefront  at  any  later  time  is  the 
envelope  of  all  such  wavelets.  These  intuitive  convictions,  sometimes 
called  Huygens'  wavefront  construction,  are  an  excellent  description  of 
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a  convolution  operation  in  which  the  initial  disturbance  is  convolved 
with  a  Huygens'  wavelet.  It  is  therefore  quite  appropriate  to  think  of 
the  point  spread  function  of  a  diffraction  system  as  the  intersection  of 
a  Huygens'  wavelet  with  the  observation  plane. 

Equation  (10)  is  therefore  an  exact  mathematical  expression  for  a 
Huygens'  wavelet  which  is  valid  right  down  to  the  initial  disturbance 
itself.  However,  for  fr  »  1, 

£2wr  i2ir  (r  -  1/4) 

h(x,y;z )  =  -i(a/r)  - -  =  (s/r)  - - 3 -  ,  (11) 

r 

it  reduces  to  the  familiar  expression  for  a  spherical  wave  with  cosine 
obliquity  factor  and  a  tr/2  phase  delay. 


General  Rayleigi  -Sommerfeld  Diffraction  Formula 

If  we  write  down  the  convolution  integral  for  the  disturbance  in 
the  observation  plane,  using  the  expression  in  Eq.  (10)  for  h(x,y;a),  we 
obtain  the  general  Rayleigh-Sommerfeld  diffraction  formula 


U{x,y\z) 


.00 


v  i2vl 

l/i*-r-dx'dy' 


(12) 


where 

i2  =  (ic-x')2  +  (y-y')2  +  a2.  (13) 

This  is  an  exact  expression  for  the  diffracted  wave  field  which  is  valid 
throughout  the  entire  space  in  which  the  diffraction  occurs--right  down 
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to  the  aperture.  No  approximations  have  been  made  in  this  scalar 
theory.  Furthermore,  the  above  equation  expresses  the  distumance  on 
the  observation  plane  as  a  superposition  of  spherical  waves  which 
corresponds  to  the  spread  functior  approach  in  image  formation. 

Geometrical  Configurations  of  the  Observation  Space 
In  order  to  insure  a  space  invariant  point  spread  function  our 
equations  have  been  restricted  to  mapping  an  optical  disturbance  from 
an  input  plane  to  an  output  plane,  where  £  has  a  parametric  relation¬ 
ship  since  its  value  determines  the  location  of  the  output  plane.  How¬ 
ever,  the  summation  of  these  Huygens'  wavelets  is  valid  over  any 
surface.  The  above  treatment  thus  gives  us  a  far  more  powerful  concep¬ 
tual  tool  than  provided  by  the  equations  themselves. 

We  will  therefore  investigate  the  properties  of  the  diffracted 
wave  field  on  two  particular  geometrical  configurations  of  the 
observation  surface. 

The  Diffracted  Wave  Field  on  a  Plane 

Equation  (12)  reduces  to  the  more  familiar  but  less  general 
form  of  the  Rayleigh-Sommerfeld  diffraction  formula  when  £  »  1, 

fr  & 

-i  K, (A’.P'jO)  4^— —  te'dg'.  (14) 

}}  t  l 

_00 

The  following  algebraic  substitutions 

1  =  *(!+«).  «  =  ^ 
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allow  us  to  rewrite  Eq.  (14)  as 
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(15) 


Note  that  we  have  imposed  no  restriction  upon  the  size  of  the  aperture 
or  the  size  of  the  observation  space.  The  only  limitation  on  the  above 
equation  is  that  the  observation  plane  must  be  many  wavelengths  from 
the  aperture. 

The  above  diffraction  formula  is  a  rather  unwieldy  integral  to 
solve  explicitly  for  most  problems  of  practical  interest.  The  Fresnel 
and  Fraunhofer  diffraction  formulas  are  obtained  by  retaining  only  the 
first  two  terms  in  the  biuomial  expansion  for  the  quantity 

l  =  g(l  +  [(£2+p2)/S2  +  (£'2+£'2)/22  -  2(££'  +PP')/S2]^. 

However,  severe  restrictions  are  then  imposed  upon  the  size  of  the 
aperture  and  the  region  over  which  the  calculations  are  valid  in  the 
observation  pla.?e.  In  order  that  we  do  not  impose  these  rest-  ictions, 
all  terms  from  the  binomial  expansion  must  be  retained.  This  can  be 
accomplished  by  rewriting  Eq.  (15)  as  the  following  Fourier  transform 
integral 

J2*z  ft  -  i|l  (x'x  +  y'y) 

U{x,y',z)  «  - —  \\%(x',9';s,y)  e  3  te'dy' ,  (16) 

iz 

-00 

where  the  complex  quantity 
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t 

* 

^(A'.P'tf.P)  =  r0(*' ,p*;0)  TT±F  ei27tf/,  (17) 

can  be  regarded  as  a  generalized  pupil  function.  TQ{&,^ ;0)  is  the  cost- 

I* 

plex  amplitude  transmittance  of  the  diffracting  aperture  (or  aperture 
function),  and  all  of  the  terras  from  the  binomial  expansion  for  the 
quantity  ( t-z ) ,  except  for  the  term  which  was  extracted  for  use  as  the 

r 

i . 

Fourier  kernel,  are  lumped  together  in  the  quantity  (f  along  with  any 
phase  variations  in  the  incident  wavefront. 

Equation  (16)  clearly  reduces  to  the  conventional  Fresnel  dif¬ 
fraction  formula  when  a  plane  wave  is  incident  upon  the  aperture  and 

A 

when  8  is  sufficiently  large  such  that  l  is  adequately  approximated  by 
retaining  only  the  first  two  terms  of  the  binomial  expansion. 

The  Diffracted  Wave  Field  on  a  Hemisphere 

Let  us  now  examine  the  diffracted  wave  field  on  a  hemisphere 
centered  upon  the  diffracting  aperture  as  illustrated  in  Fig.  2.4.  The 
position  of  an  arbitrary  observation  point  will  be  specified  by  the 
direction  cosines  a  and  6  of  its  position  vector,  and  the  radius  A  of 
the  hemisphere  upon  which  it  resides.  Note  that 

a  *  A/*,  B  =  $/*,  and  y  =  8/i\  (18) 

where 

A2  =  A2  +  02  +  82.  (19) 

The  following  algebraic  substitutions 

l  =  *(l+e);  e  =  (£-*)/*  (20) 

allow  us  to  rewrite  the  general  Rayleigh-Sommerfeld  diffraction  formula 
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We  now  have  an  exact  expression  for  the  diffracted  wave  field  on  an 
observation  hemisphere  which  is  valid  throughout  the  entire  half-space 
behind  the  plane  of  the  diffracting  aperture. 

If  we  now  require  that  £  >>  1  and  make  the  appropriate  binomial 
expansion  for  the  quantity  (£-£) ,  we  again  obtain  a  Fourier  transform 
integral 


•f.y*  *  * 


where  the  generalized  pupil  function  is  given  by 


* 

fi)  «  To{x\y';0)  jI}prTei2vW. 


Once  again,  all  of  the  terms  from  the  binomial  expansion  for  the  quan¬ 
tity  (■£-£),  except  for  the  term  which  was  extracted  for  use  as  the 
Fourier  kernel,  are  lumped  together  in  the  quantity  fc  along  with  any 
phase  variations  in  the  incident  wavefront. 
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Aberrations  of  Diffracting  Systems 

The  quantity  in  Eqs.  (17)  and  (23)  represent  phase  variations 
in  the  diffracted  wavefront  emerging  from  the  aperture.  Therefore,  ft 
can  be  interpreted  as  a  conventional  wavefront  aberration  function 
(Hopkins,  19S0)  which  is  conveniently  expressed  as  a  general  power 
series  expansion  of  the  pupil  coordinates  and  the  appropriate  field 
parameters . 

For  the  case  of  a  rotationally-symmetric  diffracting  aperture  we 
can,  without  loss  of  generality,  choose  the  observation  point  to  lie  on 
the  y- axis  (x  *=  0) .  The  wavefront  aberration  function  can  then  be 
written  as 

W  *  ^200o2  +  ^020®2  +  frlllP^  cos$ 

+  &400P**  +  ^040^i*+  Vi3ipO?COS<fc  +  ^222P2P[2 

+  ^220P2Q2  +  ^31lP3a  cos 

+  higher-order  terms,  (24) 

where  p  is  a  normalized  field  position  of  the  observation  point  and  a 
is  a  normalized  pupil  height. 

By  equating  coefficients  of  the  corresponding  terms  in  the 
appropriate  binomial  expansions  and  the  above  wavefront  aberration  func¬ 
tion,  we  obtain  expressions  for  the  aberration  coefficients  in  terms  of 
the  aperture  diameter,  the  observation  distance,  and  the  appropriate 
field  parameter.  These  aberrations,  which  are  inherently  associated  with 
the  diffraction  process,  are  precisely  the  effects  ignored  when  making 
the  usual  Fresnel  and  Fraunhofer  approximations.  Furthermore,  these 
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aberrations  have  the  same  functional  form  as  the  familiar  aberrations 
caused  by  the  refraction  or  reflection  process  in  imperfect  imaging 
systems . 

The  expressions  for  these  aberration  coefficients  are  derived 
in  Appendix  A  and  tabulated  in  Table  2.1  to  enable  easy  comparison  of 
several  different  geometrical  configurations  of  the  incident  wavefront 
and  the  observation  space. 

Consider  first  a  plane  wave  illuminating  the  diffracting  aper¬ 
ture  and  a  plane  observation  space.  We  see  from  the  first  column  of 
Table  2.1  that  all  aberrations  are  present  except  for  lateral  magnifi¬ 
cation  error  (^m)  which  is  absent  for  all  geometrical  configurations 
because  this  term  of  the  binomial  expansion  for  l  is  extracted  for  use 
as  the  kernel  of  the  Fourier  transform  integral.  It  is  clear  that  very 

A 

large  observation  distances  are  required  to  reduce  defocus  (V020)  t0  a 
negligible  value.  Also,  distortion  (V311)  imposes  severe  restraints 
upon  the  field  angle  over  which  the  diffracted  wave  field  is  accurately 
described  by  the  Fourier  transform  of  the  aperture  function.  These 
restrictions  are  the  same  as  those  usually  imposed  during  the  develop¬ 
ment  of  the  Kirchhoff  theory  and  in  most  applications  of  the  Rayleigh- 
Sommerfeld  diffraction  theory. 

The  effect  of  illuminating  the  aperture  with  a  spherical  wave 
converging  to  the  observation  plane  is  to  eliminate  defocus  (#020)  and 
all  orders  of  spherical  aberration  C^o4o) •  This  removes  the  require¬ 
ment  for  an  extremely  large  observation  distance,  but  the  Fourier 
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Table  2.1.  Tabulation  of  Expressions  for  the  Aberration  Coefficients 
for  Several  Different  Geometrical  Configurations. 

(a)  Incident  wavefront. 

(b)  Observation  space. 
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transform  of  the  aperture  function  is  still  valid  only  over  a  small 
region  about  the  optical  axis  in  the  observation  plane. 

Choosing  the  observation  space  to  be  a  hemisphere  centered  upon 
the  diffracting  aperture  eliminates  field  curvature  (^220)*  distortion 
(^3li)>  an<*  all  orders  of  piston  error  (f/200  anc*  ^4oo)- 

Hence  for  the  case  of  a  spherical  incident  wave  converging  to 
the  intersection  of  the  optical  axis  and  an  observation  hemisphere, 
only  coma  (fr131)  and  astigmatism  (V222)  are  present.  And  the  values  of 
the  aberration  coefficients  can  be  calculated  from  the  relationships 
provided  in  the  last  column  of  Table  2.1. 

Thus,  for  a  system  with  an  aperture  diameter  of  1  mm  and  an 

\ 

observation  hemisphere  with  a  radius  of  1  m,  we  have  for  X  *  0.5  wm  and 


^131  = 

1.25  x  10-4 

A 

^222  * 

-2.50  x  10-1, 

(25) 

Hence  there  is  only  A/4  of  astigmatism  at  the  edge  of  the  field  (i.e., 
90°  field  angle). 

A 

Similarly  for  an  f/10  system  ( d/r  =  0.1)  with  an  aperture  diam¬ 
eter  of  5  cm  and  a  maximum  field  angle  of  0.5°  ($max  =  0.00872),  we 
have  for  X  *  0.5  pm, 

ft131  =  0.545 

£'222  =  -0.095.  (26) 

Hence  coma  dominates  at  the  edge  of  the  field  with  a  value  of 
approximately  X/2. 
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Finally,  for  an  f/10  system  with  an  aperture  diameter  of  1  cm 
and  a  field  size  equal  to  the  size  of  the  aperture  (i.e.,  a  pair  of 
f/10  relay  lenses  1  cm  in  diameter),  we  have  for  X  =  0.5  pm, 

WUi  =  0.625 

^222  =  “0.625.  (27) 

We  find  slightly  more  than  X/2  of  coma  and  astigmatism  at  the  edge  of 
the  field. 

It  should  be  pointed  out  that  in  each  of  the  above  cases  the 

radius  of  the  Airy  disc  in  direction  cosine  space  is  approximately 

equal  to  0.001  B  .  Hence  the  off-axis  aberrations  are  of  little  more 
max 

than  academic  interest  unless  there  is  some  structure  in  the  aperture 
with  high  spatial  frequency  content  which  will  diffract  light  at  large 
angles  from  the  direction  of  the  incident  beam. 

However,  the  above  analysis  of  the  aberrations  associated  with 
the  diffraction  process  can  be  readily  applied  to  holographic  systems 
or  systems  containing  diffraction  gratings. 

For  example,  an  f/6  system  with  a  10-line-per-mm  Ronchi  ruling 
placed  in  a  40-mm  diameter  aperture  produces  the  diffraction  pattern 
shown  in  Fig.  2.5.  The  diffracted  order  at  three  different  field  posi¬ 
tions  was  photographed  through  a  microscope  with  the  following  results: 
at  6  =  0  no  aberrations  were  apparent;  at  3  =  0.04  coma  was  predominant^ 
with  a  value  of  approximately  5A;  and  at  3  =  0.10  coma  and  astigmatism 
both  have  values  of  approximately  15X.  These  values  were  estimated  by 
visual  inspection  of  the  aberrated  diffracted  orders. 


Fig.  2.5.  Diffraction  Pattern  of  a  10-Line  per  mm  Ronchi  Ruling  Placed  in 
an  f/6  Cone  of  Light  with  a  40-mm  Diameter. 

Magnified  images  of  diffracted  orders  at  various  field  positions 
indicate  that  coma  is  predominant  for  small  field  angles  with 
astigmatism  also  becoming  significant  at  larger  field  angles. 
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By  stopping  the  aperture  down  to  a  diameter  of  12  mm  (f/20 
system)  and  observing  the  diffracted  orders  at  larger  field  angles, 
astigmatism  becomes  the  predominant  aberration  as  shown  in  Fig.  2.6. 

The  sagittal  focus  lies  on  the  observation  hemisphere  of  radius  r  and 
the  medial  and  tangential  surfaces  have  smaller  radii  as  indicated. 

The  diffracted  order  at  3  8  0.020  exhibits  about  6A  of  astigmatism. 

This  order  was  observed  through  focus  with  the  microscope  and  the  mag¬ 
nified  images  are  displayed. 

In  both  of  the  above  examples  the  observed  aberrations  are  in 
good  agreement  with  those  predicted  from  the  coefficients  presented  in 
Table  2.1. 

Shift  Invariance  of  the  Diffracted  Wave  Field 
We  have  shown  that  any  departures  of  the  actual  diffracted  wave 
field  from  that  predicted  by  the  simple  Fourier  transform  of  the  aper¬ 
ture  function  take  the  form  of  conventional  wavefront  aberrations. 

If  we  neglect  these  aberrations,  Eq.  (22)  reduces  to 

U{ 0,6;*)  8  Y  ||  ^(A'^'jO)  e"i2,r(a&,+^,)d£'dp’.  (28) 

~QO 

But  this  is  merely  the  Fourier  transform  of  the  aperture  function  mul¬ 
tiplied  by  a  spherical  Huygens'  wavelet. 

t2irf 

U(a, 6;f)  8  y  Lrft-&{To&,9;0)}.  (29) 
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This  relationship  is  valid  provided  that  the  observation  space  is  a 
hemisphere  centered  on  the  diffracting  aperture  and  if  the  incident 
radiation  is  a  unit  amplitude  spherical  wave  whose  center  of  curvature 
lies  on  the  intersection  of  the  observation  hemisphere  with  the  z-axis. 
Furthermore,  if  r  is  large  compared  to  the  size  of  the  diffracting 
aperture,  the  Fourier  transform  relationship  is  accurate,  not  merely 
over  a  small  region  about  the  z-axis,  but  instead  over  the  entire 
hemisphere. 

Now  consider  the  situation  where  the  incident  radiation  strikes 
the  diffracting  aperture  at  an  angle  60  as  illustrated  in  Fig.  2.7. 

This  is  equivalent  to  introducing  a  linear  phase  variation  across  the 
aperture.  By  applying  the  shift  theorem  (Bracewell,  1965)  of  Fourier 
transform  theory  to  Eq.  (29)  we  find  that  the  complex  amplitude  distri¬ 
bution  on  the  hemisphere  is  a  function  of  (B-Bo), 

£2tt? 

y( a,3-e0;*)  =  Y  ^m*,p;0)  exp[i2v&0y]),  (30) 

where  B  is  the  direction  cosine  of  the  position  vector  of  the  observa¬ 
tion  point,  and  B0  is  the  direction  cosine  of  the  position  vector  of 
the  undiffracted  beam.  Note  that  these  direction  cosines  are  obtained 
by  merely  projecting  the  respective  points  on  the  hemisphere  back  on  to 
the  plane  of  the  aperture  and  normalizing  to  a  unit  radius.  The  com¬ 
plex  amplitude  distribution  at  an  arbitrary  point  on  the  hemisphere  can 
now  be  said  to  be  a  function  of  the  distance  of  the  observation  point 
from  the  undiffracted  beam  in  direction  cosine  space. 
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Fig.  2.7.  Geometrical  Configuration  when  the  Incident 
Beam  Strikes  the  Diffraction  Aperture  at  an 
Arbitrary  Angle. 
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As  a  specific  example,  suppose  we  have  incident  light  striking 
a  diffraction  grating  at  an  angle  60.  The  diffracted  orders  will 
strike  the  observation  hemisphere  in  a  cross  section  which  is  not  a 
great  circle  but  instead  a  latitude  slice  as  shown  in  Fig.  2.8.  Thus 
for  large  angles  of  incidence  the  various  orders  appear  to  lie  in  a 
straight  line  only  if  they  are  projected  down  onto  the  a-8  plane  in 
direction  cosine  space.  It  is  therefore  clear  that  varying  the  angle 
of  incidence  merely  shifts  the  diffracted  wave  field  in  direction 
cosine  space  without  changing  its  functional  form.  This  has  been  veri 
fied  experimentally  by  mounting  a  diffraction  grating  at  the  center  of 
a  transparent  hemisphere,  placing  graph  paper  on  the  plane  of  the  dif¬ 
fraction  grating  (a-8  plane),  and  scribing  appropriate  latitude  lines 
on  the  hemisphere  upon  which  the  diffracted  orders  fall  when  illumi¬ 
nated  with  a  small  laser  beam. 


Summary 

We  have  developed  a  very  useful  treatment  of  near-field  scalar 
diffraction  theory  that  yields  much  more  insight  than  the  conventional 
Rayleigh-Sommerfeld  theory. 

By  describing  the  diffraction  process  in  terms  of  the  direction 
cosines  of  the  propagating  light  we  have  obtained  the  extremely  power¬ 
ful  result  that  the  diffracted  wave  field  on  an  observation  hemisphere 
is  given  directly  by  the  Fourier  transform  of  the  aperture  function. 
This  allows  us  to  apply  the  well-known  techniques  of  linear  systems 
theory  that  have  proven  so  useful  in  the  area  of  image  formation. 
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Furthermore,  we  have  shown  that  any  departures  of  the  actual 
diffracted  wave  field  from  those  predicted  by  the  Fourier  transform 
relationship  take  the  form  of  conventional  aberrations  whose  behavior 
is  well  understood  in  terms  of  the  dimensions  of  the  diffraction  aper¬ 
ture,  the  radius  of  the  observation  hemisphere,  and  the  appropriate 


field  parameters. 


t 


as 


t. 


CHAPTER  3 

SURFACE  SCATTER  THEORY 

In  the  following  treatment  the  scattering  of  light  from  optical 
elements  is  considered  to  be  solely  a  surface  phenomenon.  It  is  recog¬ 
nized  that  bulk  scattering  mechanisms,  such  as  photon-phonon  interac¬ 
tions  (Bloembergen  and  Shen,  1965)  and  scattering  from  free  electrons 
(Vachaspati,  1964),  can  exist  if  the  substrate  material  is  not  perfectly 
conducting.  The  excitation  of  surface  plasmons  has  also  been  suggested 
by  several  investigators  as  contributing  to  short  wavelength  scattering 
from  polished  metal  surfaces  (Beaglehole,  1970;  Beaglehole  and  Hunderi, 
1970;  Crowell  and  Ritchie,  1970;  Elson  and  Ritchie,  1971;  and  Daude, 
Savary,  and  Robin,  1972).  However,  the  above  effects  are  believed  to 
be  small  for  roost  visible  and  infrared  radiation  scattered  from  metal 
surfaces. 

Surface  Scatter  Phenomena  as  a  Diffraction  Process 
In  Chapter  2  it  was  shown  that,  under  the  proper  circumstances, 
the  diffracted  wave  field  on  a  hemisphere  is  given  directly  by  the 
Fourier  transform  of  the  complex  amplitude  transmittance  of  the  diffrac¬ 
ting  aperture.  Usually  a  diffracting  aperture  consists  of  a  "hole"  in 
some  opaque  surface.  This  is  a  binary  amplitude  diffracting  aperture. 
Clearly,  a  continuous  amplitude  diffracting  aperture  (a  piece  of  photo¬ 
graphic  film,  for  example)  can  also  exist.  A  more  general  situation 
is  the  complex  diffracting  aperture  which  exhibits  both  amplitude  and 
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phase  variations.  It  is  these  phase  perturbations  that  lead  to  scat¬ 
tered  radiation.  Surface  scatter  phenomena  can  thus  be  described  as 
a  diffraction  process  in  which  the  pupil  function  has  random  phase 
variations  in  addition  to  any  existing  amplitude  variations.  The  dif¬ 
fraction  theory  of  the  previous  chapter  can  therefore  be  applied 
directly  to  the  problem  of  predicting  the  complex  amplitude  distribution 
on  an  observation  hemisphere  of  radius  r  resulting  from  an  incident 
beam  of  light  being  reflected  from  a  rough  surface. 

The  System  Pupil  Function 

A  simple  treatment  of  surface  scatter  theory  can  be  formulated 
by  considering  the  effect  of  the  scattering  surface  to  be  a  space- 
dependent  modifier,  or  random  component,  of  the  effective  pupil  func¬ 
tion  of  the  system.  The  disturbance  emerging  from  the  scattering  sur¬ 
face  is  then  given  by 

o)  =  aoPC&.p;  o)  =  ajp^x^^p^x^-yO).  (3i) 

Here  the  pupil  function  of  the  system  producing  the  incident 
beam  is  given  by 

i2vWL(x,y;0) 

0}  *  aL(x,y;0)  e  ,  (32) 

where  a ^  describes  the  amplitude  variations  across  the  exit  pupil  of 
the  system,  and  WL  describes  any  phase  variations  or  aberrations  in 
the  wavefront  of  the  incident  beam. 

The  random  component  of  the  pupil  function  due  to  the  scatter¬ 
ing  surface  similarly  has  an  amplitude  and  phase  component 


Pff(x>9‘>  o) 


M&,9;o j 


(33) 


e 

Here  i?  is  the  reflectance  of  the  scattering  surface  and  Wg  is  the  sur¬ 
face  height.  Note  that  the  phase  variations  on  the  reflected  wavefront 
are  twice  as  large  as  the  actual  variations  on  the  reflecting  surface. 

A 

Figure  3.1  illustrates  the  surface  heigl  t  variations,  V%,  as  a 
function  of  distance  along  the  surface.  This  surface  profile  has 
associated  with  it  an  autocovariance  function  and  a  surface  height 
distribution  function  as  indicated. 


Surface 

height  Surface 

distribution  height 


Fig.  3.1.  Illustration  of  Surface  Height  Variations  and  Associated 
Statistical  Parameters. 
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Intensity  Distribution  of  a  Scattering  System 

The  results  of  the  previous  chapter,  [Eq.  (29)],  indicate  that 
the  complex  amplitude  distribution  on  the  observation  hemisphere  is 
given  directly  by  the  Fourier  transform  of  U0{&,9 ;0) 

i2-nS> 

0(<*,e;fO  *  Y  “ — ~  0{Uo@,9i 0)}.  (34) 

ir  * 

The  total  reflected  flux  $  is  obtained  by  applying  Rayleigh's 
theorem  from  Fourier  transform  theory 

00  00 

*  *  |[  |2<iad$  *  a*  ||f(A',P';0) \2<&'dy'. 

(35) 

Noting  that  dta  »  dctd&/y,  the  radiant  intensity  of  the  scattering 
system  can  be  written  as 

-T(o.8)  *  gjj’  s  “  |tf(o,S;r)  |2  *  ao2y\^{p(x,yi0)}\z. 

(36) 

Utilizing  the  autocorrelation  theorem  of  Fourier  transform  theory,  this 
is  equivalent  to 

|j  P@\9';0)p*(x'-£,y'-y;0)dx'dy^.  (37) 

For  the  special  case  of  a  Lambertian  surface,  the  autocorrela¬ 
tion  function  approaches  a  Dirac  6-function;  hence,  its  Fourier  trans¬ 
form  is  constant  and  we  obtain 

IL(ot0)  »  y K,  K  *  constant  (38) 

which  is  consistent  with  Lambert's  cosine  law. 


38 


The  System  Spread  Function 

Following  the  standard  procedure  used  in  image  evaluation,  the 
effective  transfer  function  of  the  scattering  system  is  defined  as  the 
normalized  autocorrelation  of  the  pupil  function 


oo 

||  p(£'»i/';< 


0)p*$'-x,y'-y;0)dx'dy' 


* 

|p(A’,p';0)|2 


(39) 


The  effective  spread  function  of  the  scattering  system  is  now 
defined  in  the  usual  way  as  the  Fourier  transform  of  the  transfer 
function 

/  °° 

0x||  P(*\p^0)p*(a'-*,p'-p;0)cfc'dp'j 

j(a,8)  =  -  — — _ _ _ l 

00 

II  |p(£',P';0)|2c&'d&'  (40) 


-00 


Direct  substitution  from  Eq.  (35)  and  Eq.  (38)  results  in  the  following 
expression  for  the  effective  spread  function  in  t/rms  of  the  radiant 
intensity  of  the  scattering  system 


/ 

\ 


<2(a,B) 


J(«. 6) 

$(l-a2-g2)* 


(41) 


Scattering  from  Optical  Surfaces 
Let  the  height  variations  of  a  given  illuminated  area  on  a 
scattering  surface  be  a  two-dimensional  sample  function  A 
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A  A  ^ 

random  process,  W^{x,y;0),  is  made  up  of  an  ensemble  of  such  functions 
as  shown  in  Fig.  3.2.  For  fixed  spatial  coordinates,  j,pi;0)  is  a 
random  variable.  And,  for  a  specific  sample  function  with  fixed 
spatial  coordinates,  &L,(£1,p1;0)  is  a  single  number. 


Wi(x,y) 


>,.*» ! 

jW2(x>yl  j 

U  i 


Cx2,y2)> 


;Cx2  .y23 , 


i(xi.yi)  | 

i* . .  i 

X(x,y)  . 


I(X2.y2)i 


](xj,y2)  j 

i — f 


'"'(xj.yi) 


W(x^,y2) 


Fig.  3.2.  An  Ensemble  of  Two-dimensional  Sample  Functions  #£(&,*/) 

Representing  Surface  Height  Variations  Constitutes  a  Random 
Process  fr(2,p). 

Two  random  variables,  and  W(&2,y 2)  with  fixed 

spatial  coordinates  are  also  shown. 
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The  Statistical  Properties  of  an  Optical  Surface 

We  are  primarily  interested  in  the  scattering  behavior  of  opti¬ 
cal  surfaces.  The  following  assumptions  are  made  concerning  the  sta¬ 
tistical  properties  of  an  optical  surface  prepared  by  conventional 
fabrication  techniques  on  ordinary  optical  mateiials: 

1.  The  reflectance  R  is  constant  over  the  entire  surface.  This 
assumption  is  not  essential  but  it  is  reasonable  and  furthermore  it  has 
been  shown  by  Shack  (1967)  that  phase  fluctuations  will  dominate  over 
amplitude  fluctuations  in  their  combined  effect  on  the  spread  function. 

2.  0)  is  a  single-valued  Gaussian  random  process. 

3.  0)  is  at  least  locally  stationary  in  the  statistical 
sense  (i.e.,  surface  is  homogeneous  and  isotropic). 

4.  The  random  variables  ,pi  ;0)  and  ^>(^2»P 2!°)»  produced  by 

any  two  fixed  pairs  of  spatial  coordinates,  are  jointly  normal. 

5.  0)  is  weakly  ergodic  (i.e.,  the  mean  and  autocorrela¬ 

tion  function  determined  by  space  averages  using  a  single  sample  func¬ 
tion  0)  are  the  same  as  those  determined  by  ensemble  averages). 
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The  Transfer  Function  of  a  Scattering  Surface 

Substituting  Eq.  (31)  into  Eq.  (39)  we  obtain 


00 

II 


A  A 


i2v[WLi-WL2]  i4*[Wni-WR2] 
aLlaL2  e  e  dx'dy' 


Jf(xty;0)  = 


||  \aL1\z  dx'dy' 


(42) 


where 


aLl 

Ujr  (ic*  *  ;0) 

aL2 

= 

o  (a1-®,  y'-y; 

0) 

A 

hi 

s 

h2 

= 

Wl(x'-3:,  y'-y; 

0) 

hi 

= 

WR(x',y",  0) 

h2 

= 

^(x’-x,  y’-y; 

0). 

(43) 


The  above  expression  for  the  transfer  function  contains  the  ran¬ 
dom  variables  and  therefore,  taking  the  expected  value  we  have 


w 

11“ 


•LlaL2e 


E{J tr(x,y;0)}  = 


i2v[WLi-WL2]  (  t4Tr[Vi?rfc'i?2]) 

9  F<e  >dx’d£’ 


(44) 


||  \aL\2dx'dQ' 
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that  (see  Papoulis,  1965,  p.  226) 


£<e  >  =  e" 


4Tr(ni-n2)e-87r2  (oi2-2Ci2+o22) 


(48) 


where 


C12  =  ^{fe-ni)  C^i?2-n2)} 


(49) 


is  the  covariance  function  of  the  random  variables  and  But  W ^ 

A 

and  are  identical  functions  merely  displaced  from  one  another;  hence 


01 


a2 


ni 


n2 


(50) 


and 


C12  -  cfy(x,y)  =  autocovariance  of  W^. 


(51) 


The  equivalent  transfer  function  of  the  scattering  surface  is  thus  given 
by 


-  (4*0#)  1 


j*^(x,p;0)  *  e 


1  - 


Cp  (*.$) 

—ZTt — 
aW 


(52) 


where  o^2  is  the  variance  of  the  surface  height  distribution  function 
and  C#(£»p)  is  the  two-dimensional  autocovariance  function  of  the 
surface. 

Considerable  insight  into  the  scattering  process  can  now  be 
obtained  by  considering  the  nature  of  this  transfer  function.  The  auto¬ 
covariance  function  approaches  the  value  o#2  as  the  displacement 
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approaches  zero.  The  equivalent  transfer  function  thus  approaches  unity 
as  expected.  As  the  displacement  approaches  infinity  the  autocovariance 
function  approaches  zero  and  the  equivalent  transfer  function 
approaches  a  plateau  of  height  exp[-(2ira£0*]  • 

The  equivalent  transfer  function  of  the  scattering  surface  can 
thus  be  regarded  as  the  sum  of  a  constant  component  and  a  bell-shaped 
component  as  shown  in  Fig.  3.3(a).  Equation  (52)  can  therefore  be 


rewritten  as 


Jt£Cx,£;0)  «  A*B-Q{&,P; 0), 


where 


A  *  e 


■  (4irOy): 


B  -  1  -  e 


-(4no£r): 


(4TTOy) 


.  o-jM) 

— _ 


Q(x,p;  0) 


(4<TOy): 


The  Spread  Function  of  a  Scattering  Surface 

The  significance  of  this  interpretation  of  the  equivalent  trans¬ 
fer  function  of  the  scattering  surface  is  dramatically  shown  by  the 
inferred  properties  of  the  corresponding  spread  function.  Since  the 
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transfer  'function  is  the  sum  of  two  separate  components,  the  equivalent 
spread  function  of  the  scattering  surface  is  the  sum  of  the  inverse 
Fourier  transforms  of  the  two  component  functions. 


where 


J(a ,$;*) 


5(a,g;r) 


*  =  AS  (a, B;r)  +  S(a,M) 


(57) 


The  constant  component  transforms  into  a  delta  function,  and  the  bell¬ 
shaped  component  transforms  into  a  bell-shaped  scattering  function  as 
shown  in  Fig.  3.3(b).  Hence  the  scattering  surface  reflects  an  incident 
beam  of  light  as  a  specularly-reflected  beam  of  diminished  intensity 
surrounded  by  a  halo  of  scattered  light.  Furthermore,  the  relative 
power  distribution  between  the  specular  component  and  the  scattered 
component  of  the  effective  spread  function  are  given  by  the  quantities 
A  and  B  respectively. 

Note  that  as  more  and  more  light  is  scattered,  energy  is 
transferred  from  the  specular  component  of  the  spread  function  into 
the  scattered  component  of  the  spread  function.  For  a  perfectly 
Lambertian  reflector  the  specular  component  disappears  completely 
from  the  spread  function. 


A  Shi ft- invariant  Scattering  Function 

In  general,  the  scattered  light  distribution  on  an  observa¬ 
tion  hemisphere  will  change  with  the  angle  of  the  incident  light  just 
as  the  point  spread  function  of  an  optical  imaging  system  will,  in 
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general,  vary  with  the  field  position  of  the  point  source.  However, 
the  analysis  of  imaging  systems  is  greatly  simplified  by  assuming  an 
isoplanatic  system  in  which  the  point  spread  function  does  not  change 
with  field  position  (and  this  is  a  reasonable  assumption  for  many 
practical  imaging  systems).  Similarly,  the  analysis  of  light  scatter¬ 
ing  systems  will  be  greatly  simplified  if  they  can  be  shown  to  be 
shift- invariant  (i.e.,  if  the  shape  of  the  scattering  function  does 
not  change  with  the  angle  of  incidence) . 

From  the  discussion  in  the  previous  section  it  is  clear  that 
the  scattered  light  distribution  on  an  observation  hemisphere  will 
appear  to  consist  of  the  sum  of  two  components,  a  core  which  is  the 
delta  function  convolved  with  the  spread  function  of  the  optical 
system  producing  the  incident  beam,  and  a  scattering  function  which 
is  the  bell-shaped  halo  convolved  with  the  spread  function  of  the 
optical  system. 

In  Fig.  3.4  we  have  merely  replaced  the  diffracting  aperture 
of  Fig.  2.7,  page  31,  with  a  scattering  surface  and  the  geometry  of 
the  measurements  has  been  folded  about  the  reflecting  plane.  Hence, 
we  have  the  incident  beam  striking  the  scattered  surface  at  some  angle 
of  incidence,  a  specularly-reflected  beam  striking  the  observation 
hemisphere,  and  the  scattered  light  distribution  being  sampled  at  an 
.arbitrary  point  with  direction  cosine  coordinates  a  and  6.  The  scat¬ 
tered  light  distribution  on  the  hemisphere  will,  in  general,  change 
shape  drastically  with  angle  of  incidence--becoming  quite  skewed  and 
asymmetrical  at  large  angles  of  incidence.  However,  our  theory 
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Fig.  3.4.  Geometrical  Relationship  Between  Incident  Beam,  Scatter¬ 
ing  Surface,  the  Measured  Scattered  Light  Distribution, 
and  the  Resulting  Spread  Function. 
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indicates  that  the  data  collected  on  the  hemisphere  should  be  plotted 
as  a  function  of  the  direction  cosines  of  the  position  vector  of  the 
observation  point.  For  certain  surfaces  with  well-behaved  statistics, 

& 

this  new  scattering  function  will  not  change  shape  but  will  merely  be 
shifted  in  direction  cosine  space  with  changes  in  angle  of  incidence. 
The  four- dimensional  Bidirectional  Reflectance  Distribution  Function 

J 

(BRDF),  which  is  the  basic  quantity  that  completely  characterizes  the 
scattering  properties  of  a  surface  (see  Appendix  B),  will  therefore 
degenerate  into  a  single  two-dimensional  spread  function.  This  is  a 
rather  significant  development  which  has  profound  implications  regard¬ 
ing  the  quantity  of  data  required  to  completely  characterize  a  scatter¬ 
ing  surface.  However,  it  remains  to  be  experimentally  verified  that 
scattering  surfaces  of  practical  interest  obey  these  predictions. 

The  Inverse  Scattering  Problem 

The  problem  of  determining  surface  characteristics  from  scat¬ 
tered  light  measurements  is  frequently  referred  to  as  the  inverse 
scattering  problem.  A  general  treatment  of  electromagnetic  inverse 
scattering  has  been  discussed  by  Bojarski  (1971).  Several  attempts 
have  been  made  (Daniels,  1961;  Fung  and  Moore,  1964;  and  Barrick, 

1965)  to  determine  properties  of  the  lunar  surface  by  applying 
inverse  scattering  techniques  to  radar  returns  from  the  moon.  Bennett 
and  Bennett  (1967)  were  able  to  obtain  the  rms  roughness  and  autocovar¬ 
iance  length  of  the  surface  structure  of  mirrors  by  assuming  a  Gaussian 
shape  for  both  the  surface  height  distribution  function  and  the  surface 
autocovariance  function.  More  recently  Scheele  (1973)  met  with  little 
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success  in  attempting  to  ascertain  under  what  conditions  the  exact 
autocovariance  function  can  be  obtained  from  scattered  light  data. 
The  treatment  presented  in  this  chapter  describes  surface 


scatter  phenomena  as  a  linear,  shift-invariant  process  which  is 
completely  characterized  by  the  effective  transfer  function  of  the 
scattering  surface.  Furthermore,  for  a  large  class  of  well-behaved 
surfaces  this  transfer  function  is  described  only  in  terms  cf  the  rms 
surface  roughness  and  the  surface  autocovariance  function  and  hence 
provides  an  elegant  solution  to  the  inverse  scattering  problem. 

The  surface  autocovariance  function  is  thus  obtained  from 
scattered  light  data  by  rewriting  Eq.  (52)  as 

I 

'  (&)*  *" (*«(!'  !))“*  (58) 

where 

*&'*)  |*.*/>  (59) 

y  *  j/A  . 

The  rms  surface  roughness  is  given  by 


where  B  (the  total  integrated  scatter  expressed  as  a  fraction  of  the 
total  reflected  light)  is  obtained  by  applying  numerical  integration 
techniques  upon  the  measured  data  describing  the  scattered  light  dis¬ 
tribution.  Note  that  although  we  are  limited  by  a  scalar  theoTy,  we 


have  made  no  explicit  approximations  regarding  the  size  of  the  surface 
variations.  If  we  make  a  smooth  surface  approximation  (<fy«A),  then 
the  total  integrated  scatter  is  small  (£«1)  and  we  obtain  the  usual 
expression  for  the  rms  surface  roughness 


aW 


(61) 


Wavelength  Dependence  of  the  Scattering  Function 

A  successful  theoretical  model  of  surface  scatter  phenomena 
must  provide  a  method  of  determining  the  wavelength  dependence  of  the  ' 
scattered  light  characteristics.  This  would  allow  one  to  infer  the 
scatteiing  behavior  of  a  surface  for  any  desired  wavelength  from  a 
limited  amount  of  data  obtained  at  a  given  wavelength.  Since  the  scat¬ 
tering  mechanism  is  a  diffraction  process,  it  is  clear  that  light  of  a 
particular  wavelength  scattered  in  a  given  direction  corresponds  to 
surface  structure  of  a  given  spatial  frequency.  This  spatial  frequency 
component  of  the  surface  structure  will  scatter  light  of  some  other 
wavelength  into  a  different  direction.  An  analytical  expression  de¬ 
scribing  the  wavelength  dependence  of  the  scattering  function  must 
therefore  involve  a  change  in  the  scatter  angle  with  wavelength  as 
well  as  the  expected  change  in  scattered  intensity.  Thus,  if  we 
wish  to  determine  how  the  relative  intensity  at  a  fixed  scatter  angle 
varies  with  wavelength,  the  surface  characteristics  corresponding  to 
the  appropriate  spatial  frequency  components  of  the  surface  structure 
must  be  known  for  each  wavelength  of  interest. 
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The  above  discussion  perhaps  explains  the  failure  of  various 
experimental  investigators  to  agree  in  their  attempts  to  establish  the 
wavelength  scaling  behavior  from  empirical  observations  of  scattered 
light  data.  Shack  and  DeBell  (1974)  made  scattered  light  measurements 
on  mirror  surfaces  at  two  visible  wavelengths.  Their  data  indicated 
a  1/A2  wavelength  scaling  law.  Leinert  and  Kltlpelberg  (1974)  also  made 
mirror  scatter  measurements  at  two  wavelengths  in  the  visible.  They 
found  a  1/A  wavelength  dependence.  Perkin-Elmer  (1975)  mirror  scatter 
measurements  were  made  at  a  visible  wavelength  and  at  A=10ym.  This 
data  showed  a  1/A^  wavelength  scaling  behavior. 

The  transfer  function  characterization  of  scattered  surfaces 
developed  in  the  previous  section  offers  a  simple  means  of  determining 
this  wavelength  dependence.  Once  the  surface  characteristics  are  known 
[whether  from  direct  measurement  or  calculated  from  scattered  light  data 
by  means  of  Eq.  (58)  and  Eq.  (60)],  the  same  theory  provides  an  equally 
simple  method  of  predicting  the  scattering  function  at  any  desired  wave¬ 
length. 

Since  this  technique  involves  numerical  computations  on  sampled 
data,  an  analytical  expression  for  the  wavelength  scaling  law  is  not 
required.  However,  in  order  to  gain  more  insight  into  the  nature  of 
surface  scatter  phenomena,  we  will  proceed  to  derive  the  wavelength 
scaling  law  for  the  special  case  of  a  normally  incident  beam  upon  a 
smooth  surface  (cty<< A).  Under  this  condition  the  surface  transfer  func¬ 
tion  expressed  in  terras  of  real  pupil  coordinates  x  and  y  is  given  by 
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where 


RR(x,y)  *  A  +  BQ(x,y ) 


>1  -  1 


-p) 

py 


Q(*>y) 


Cw(x,y) 


From  Eq.  (57),  we  see  that  the  scattering  function  for  a  particular 
wavelength  is  given  by 


S(a,B;X)  -  A*,A0)}. 


Applying  the  similarity  theorem  of  Fourier  transform  theory  we  obtain 


where 


SC«.M)  =  (t)2  iWf  •  l) 


^(a.8)  «  rf-1  Cw(x,y)  . 


If  we  r.ov*  scale  the  wavelength  by  a  factor  a,  we  obtain 

SM-.aX)  .  (il)2  ^  .  &)  •  Jr  ^  •  f  i»).  (69) 

Therefore,  the  appropriate  wavelength  scaling  law  for  smooth  surfaces 
is  given  by 

S(a,8;aA)  *  s(£  ,  |  ;a).  (70) 

Note  that,  in  addition  to  the  1/a4  change  in  magnitude,  the  width  of  the 
scattering  function  in  direction  cosine  space  is  scaled  by  the  factor  a. 
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CHAPTER  4 


SURFACE  SCATTER  MEASUREMENTS 

Apparatus 

An  instrument  has  been  designed  and  built  at  the  Optical 
Sciences  Center  for  making  scattered  light  measurements  on  a  hemisphere 
as  described  in  the  previous  chapter.  A  schematic  diagram  of  this 
apparatus  is  shown  in  Fig.  4.1.  The  incident  light  passes  through  a 
chopper  so  that  synchronous  detection  with  a  PAR  lock-in  amplifier  can 
be  made.  The  mechanical  apparatus  shown  in  Fig.  4.2  is  located  in  a 
small  photometric  darkroom  in  which  the  experiment  is  conducted.  A 
movable  arm  with  folding  mirrors  can  be  positioned  to  direct  the  inci¬ 
dent  beam  onto  the  sample  at  any  desired  angle.  A  lens  positioned  on 
this  arm  focuses  the  incident  radiation  onto  the  hemisphere  mapped  out 
by  the  detector;  hence,  the  geometrical  configuration  is  consistent 
with  that  illustrated  in  Fig.  3.4  of  the  previous  chapter.  Two  sepa¬ 
rate  driving  mechanisms  allow  us  to  measure  the  scattered  light  distri¬ 
bution  over  the  entire  hemisphere  bounded  by  the  plane  of  the  sample. 

The  Light  Source 

The  light  source  employed  is  a  Spectra- Physics  Model  165  Argon 
Ion  Laser.  The  laser  is  operated  with  a  light-regulated,  single¬ 
frequency  output  which  assures  intensity  regulation  to  within 
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Fig.  4.1.  Schematic  Diagram  of  Scatter  Measurement  Apparatus. 
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ig.  4.2.  Apparatus  for  Measuring  Scattered  Light  Distribution  from 
Optical  Surfaces. 

(1)  Precision  rotary  table,  (2)  worm  gear  drive  for  arm 
supporting  detector,  (3)  photomultiplier  tube  and  fiber¬ 
optic  probe,  (4)  sample  holder,  (S)  movable  arm  with  fold¬ 
ing  mirrors  and  lens  for  directing  and  focusing  incident 
beam  (6)  P.A.R.  lock-in  amplifier,  (7)  high  voltage  power 
supply  for  PMT. 


one  percent.  The  measurements  were  made  with  approximately  20  mw  of 
power  in  the  incident  beam  at  wavelengths  of  0.5145  ym  and  0.4579  ym. 

An  He-Ne  laser  was  also  available  for  making  measurements  at  a  wave¬ 
length  of  0.6328  ym. 

The  Detector  Unit 

The  detector  is  a  Phillips  one-inch,  end-on  photomultiplier 
tube  (PMT)  having  an  S-20  photocathode.  Light  reaches  the  photomulti¬ 
plier  by  way  of  a  rigid  fiber-optic  probe.  Such  a  probe  offers  several 
distinct  advantages  in  light  sampling.  In  addition  to  allowing 
increased  angular  resolution  throughout  the  sampling  space,  and 
enabling  us  to  sample  within  one  degree  of  the  incident  or  specu¬ 
larly  reflected  beams,  it  provides  the  ability  to  control  the  field 
of  view  of  the  detector  for  the  purpose  of  stray  light  rejection. 

The  original  configuration  consisted  of  a  rigid  fiber-optic 
bundle  bent  such  that  one  end  was  pointed  toward  the  illuminated 
spot  on  the  sample.  This  end  of  the  bundle  thus  acted  as  the  col¬ 
lecting  aperture  for  the  detection  system.  The  other  end  of  the 
fiber-optic  bundle  protruded  into  the  photomultiplier  tube  housing 
followed  by  a  series  of  baffles  to  limit  the  field  of  view  of  the 
detector  as  shown  in  Fig.  4.3(a).  This  resulted  in  a  detector 
response  which  had  a  Gaussian  dependence  upon  field  angle. 

However,  by  introducing  a  small  collecting  lens  and  a  field 
stop  in  front  of  the  fiber-optic  bundle  as  shown  in  Fig.  4.3(b),  the 
baffles  can  be  eliminated  and  a  well-defined  field  of  view  of  any 
desired  size  can  be  obtained  by  properly  choosing  the  size  of  the 


Fig.  4.3.  Illustration  of  Detector  Probe  Unit. 

(a)  Previous  Configuration. 

(b)  New  Configuration. 
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field  stop.  This  is  more  clearly  illustrated  in  Fig.  4.4.  A  coated 
doublet  with  a  10-mm  focal  length  was  edged  down  and  mounted  in  a 
black  anodized  brass  tube  3-mm  in  diameter.  A  field  stop  allowing 
a  5-degree  field  of  view  was  fabricated  and  inserted  into  the  tube 
at  the  rear  focal  plane  of  the  lens.  This  assembly  was  then  posi¬ 
tioned  onto  the  end  of  the  fiber-optic  bundle.  The  detector  response 
from  a  small  (point  source)  light  source  was  then  recorded  as  a 
function  of  fioid  angle  for  both  detector  probe  configurations. 

The  results  are  displayed  for  comparison  in  Fig.  4.5.  Both  the  flat 
response  and  the  sharp  cutoff  obtained  with  the  modified  unit  are 
highly  desirable  features.  The  flat  response  promises  to  eliminate 
signal  variations  due  to  slight  misalignment  or  wobble  ir.  the  mechani¬ 
cal  instrument  while  scanning  over  the  hemisphere.  The  ability 
to  keep  the  field  of  view  small  with  a  very  sharp  cutoff  is  essential 
for  stray  light  rejection. 

The  scattered  light  flux  from  a  polished  surface  varies  by 
several  orders  of  magnitude  over  the  angular  range  to  be  measured. 
Hence  the  linearity  of  the  PMT  was  measured  using  a  calibrated  neu¬ 
tral  density  wedge  and  several  known  neutral  density  filters  to  vary 
the  incident  flux.  The  resulting  linearity  curve  is  shown  in  Fig. 

4.6  and  indicates  a  deviation  of  less  than  1%  over  a  range  of  five 
orders  of  magnitude  of  the  incident  flux. 
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Field  angle  (degrees) 

Fig.  4.5.  Detector  Response  as  a  Function  of  Field  Angle  for  Detector 
Configurations  Shown  in  Fig.  4.3. 

(a)  Detector  configuration  shown  in  Fig.  4.3(a). 

(b)  Detector  configuration  shown  in  Fig.  4.3(b). 
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The  Scanning  Mechanism 

The  mechanical  apparatus  for  measuring  the  angular  distribu¬ 
tion  of  light  scattered  from  optical  surfaces  was  shown  in  Fig.  4.2. 

The  detector  probe  unit  is  mounted  on  a  rigid  arm  that  can  be  rotated 
in  either  of  two  orthogonal  directions.  Rotation  about  a  vertical 
axis  is  accomplished  by  means  of  a  massive  precision  rotary  table. 

The  rigid  arm  is  attached  to  the  rotary  table  by  means  of  a  worm  gear 
arrangement  that  allows  rotation  about  the  horizontal  axis.  These 
two  separate  driving  mechanisms  allow  us  to  measure  the  scattered 
light  distribution  over  the  entire  hemisphere  bounded  by  the  plane  of 
the  sample.  However,  in  order  to  limit  the  quantity  of  data  to  be 
collected,  the  scattered  radiation  field  is  sampled  in  two  principal 
directions.  These  include  the  entire  plane  of  incidence  and  a  plane 
perpendicular  to  both  the  plane  of  incidence  and  the  plane  of  the 
sample  which  also  passes  through  the  intersection  of  the  specular  beam 
with  the  observation  hemisphere  (see  Fig.  4.7).  This  particular  sam¬ 
pling  procedure  was  chosen  because  each  sampling  direction  then  involves 
one  fixed  coordinate  in  direction  cosine  space.  Furthermore,  the 
apparatus  was  designed  such  that  each  of  the  two  independent  drive 
mechanisms  corresponds  directly  to  a  given  coordinate  in  direction 
cosine  space.  Hence  for  a  given  observation  point  determined  by  the 
angles  0  and  (j>  displayed  on  the  apparatus,  the  corresponding  coordinates 
in  direction  cosine  space  are  given  by 

a  =  cos0  sin<f> 

8  =  sin0.  (71) 
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The  Incident  Beam 


A  large  movable  arm  with  appropriate  folding  mirrors  is  used 
to  direct  the  incident  beam  onto  the  sample  at  any  desired  angle  of 
incidence.  The  original  configuration  included  a  single  lens  posi¬ 
tioned  on  this  arm  to  focus  the  incident  radiation  onto  the  hemisphere 
mapped  out  by  the  detector  as  shown  in  Fig.  4.8.  This  configuration 
allowed  scattered  light  from  the  folding  mirror  M3  to  reach  the  observa¬ 
tion  hemisphere  after  being  reflected  from  the  test  sample.  The  scat¬ 
tered  light  level  from  M3  often  exceeded  that  from  the  test  sample, 
thus  rendering  the  scattered  light  data  at  small  observation  angles 
completely  useless  (Shack  and  DeBell,  1974). 

An  improved  configuration  is  shown  in  Fig.  4.9.  Immediately 
following  the  last  folding  mirror  the  beam  is  focused  onto  a  pinhole 
which  acts  as  a  spatial  filter.  This  spatial  filter  assembly  elimi¬ 
nates  from  the  beam  incident  upon  the  szmple  any  light  scattered  from 
the  folding  mirrors  as  well  as  any  diffraction  effects  from  the  chopper 
blade.  Lens  L2  then  forms  an  image  of  the  pinhole  upon  the  hemisphere 
mapped  out  by  the  collecting  aperture  of  the  scanning  fiber-optic 
probe.  The  dramatic  improvement  obtained  with  this  configuration  is 
illustrated  by  ;omparing  Fig.  4.10(a)  with  Fig.  4.10(b).  These  photo¬ 
graphs  were  obtained  by  placing  a  piece  of  photographic  film  in  the 
observation  space  at  position  B  in  Fig.  4.9  along  with  a  small  obstruc¬ 
tion  to  block  the  specular  beam.  The  six  bright  spots  in  Fig.  4.10(b) 
were  the  weld  marks  on  the  back  side  of  the  pinhole  which  were  illumi¬ 
nated  by  the  light  reflected  from  L2  then  imaged  by  L2  onto  the  observa¬ 
tion  space. 
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Fig.  4.9.  Improved  Configuration  of  Beam-forming  Optics. 
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Fig.  4.10.  Photograph  of  Stray  Light  Reaching  the  Observation  Space 
for: 


(a) 

(b) 

(c) 


Original 

Improved 

Improved 


configuration  of  beam-forming  optics, 
configuration  of  beam-forming  optics, 
configuration  with  black  pinhole. 
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The  acquisition  of  a  pinhole  in  a  black  substrate  with  no  visi¬ 
ble  weld  joints  resulted  in  the  photograph  displayed  in  Fig.  4.10(c). 
Based  upon  measurements  made  with  the  instrument  with  no  sample  in 
place,  the  stray  light  in  the  observation  space  has  thus  been  reduced 
by  almost  two  orders  of  magnitude. 

The  aberrations  associated  with  the  scattered  light  distribu¬ 
tion  on  the  observation  hemisphere  were  considered  in  determining  the 
geometrical  parameters  of  the  incident  beam  and  the  detector  scanning 
mechanism.  The  aberration  coefficients  presented  in  the  last  column  of 
Table  2.1  yield  the  following  values  for  coma  and  astigmatism  at  the 
edge  of  the  field  (0  =  90°)  when  a  spot  4  mm  in  diameter  is  illuminated 
on  the  sample  and  the  observation  distance  is  250  mm, 

fr131  =  0.128 

W222  =  16.  (72) 

Although  a  substantial  amount  of  astigmatism  exists  at  the  edge  of  the 
field,  our  tolerance  is  quite  loose  as  we  have  a  3  mm  diameter  collec¬ 
ting  aperture  on  the  fiber-optic  probe. 

Experimental  Procedures 

The  goal  of  the  research  reported  here  was  to  examine  the 
scattering  properties  of  samples  whose  surface  characteristics  span 
those  typically  produced  with  optical  fabrication  techniques. 
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Sample  Preparation 

Surface  preparation  techniques  used  to  produce  a  set  of  samples 
are  outlined  in  Table  4.1.  All  samples  were  finished  to  be  nominally 
flat. 

The  prepared  samples  were  cleaned  prior  to  coating  with  alumi¬ 
num.  Cleaning  consisted  of  careful  washing  with  Liquinox,  a  mild 
detergent,  under  very  warm,  filtered  tap  water.  Samples  were  then 
mounted  in  a  sample  holder  while  held  in  distilled  water.  Once  in  a 
holder,  samples  were  moved  to  an  ultrasonic  cleaner  filled  with  dis¬ 
tilled  water  for  rinsing.  Once  rinsed  the  samples  were  set  to  dry  in  a 
dust-free  atmosphere.  Dry  samples  were  removed  from  the  sample  holders 
and  placed  in  individual  boxes  ard  supported  by  the  edge  of  their  back¬ 
side.  Mott  (1971)  used  a  similar  cleaning  technique,  which  he 
describes  more  completely. 

Cleaned  dry  samples  were  then  placed  in  a  high  vacuum  chamber 
and  coated  to  near  opacity  with  pure  aluminum.  The  coating  technique 
varied  from  standard  procedures  only  in  that  excessive  care  was  taken 
to  allow  the  chamber  to  reach  a  pressure  below  2  *  10" 6  torr  prior  to 
coating.  The  samples  were  allowed  to  cool  to  room  temperature  prior  to 
removal  from  the  chamber.  Each  coating  run  contained  ten  different 
samples.  Once  coated,  samples  were  returned  to  their  individual  storage 
boxes.  After  all  samples  were  coated,  the  best  samples  of  each  type 
were  selected  for  measurement.  This  selection  was  made  on  the  basis  of 
individual  inspection  of  each  sample  while  held  under  a  microscope 
illuminator  in  an  otherwise  dark  room.  Samples  that  had  coating  non¬ 
uniformities,  sleeks  or  pinholes  were  rejected  as  were  those  with 


Table  4.1.  Surface  Preparation  Techniques 
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waterspots,  large  scratches,  or  otherwise  questionable  appearance. 

Prior  to  each  set  of  scatter  measurements,  samples  were  again  individu¬ 
ally  inspected  for  flaws.  Dust  was  removed  using  a  commercially  avail¬ 
able  pressurized  air  can.  After  each  sample  was  measured  for  scatter, 
it  was  returned  to  its  individual  box. 

System  Alignment 

Before  any  meaningful  scatter  measurements  could  be  made  it  was 
necessary  to  systematically  align  the  entire  system. 

The  incident  laser  beam  was  first  adjusted  to  lie  in  a  horizon¬ 
tal  plane.  Then  the  mechanical  apparatus  was  positioned  such  that  the 
axis  of  rotation  of  the  movable  arm  supporting  the  folding  mirrors  (see 
Fig.  4.9)  was  colinear  with  the  incident  beam.  This  was  accomplished 
by  means  of  four  massive  leveling  screws  at  the  base  of  the  stand  sup¬ 
porting  the  entire  apparatus. 

With  mirror  Ml  removed  the  laser  beam  passed  through  the  small 
hole  in  the  center  of  the  bearings  upon  which  the  movable  arm  rotates. 
This  assured  that  the  incident  beam  was  indeed  colinear  with  the  axis 
of  rotation  and  furthermore  allowed  the  sample  holder,  which  was 
mounted  on  a  shaft  passing  througli  the  axis  of  the  precision  rotary 
table,  to  be  accurately  positioned  with  the  center  of  the  scattering 
surface  lying  at  the  intersection  of  these  two  perpendicular  axes.  The 
precision  rotary  table  was  then  positioned  such  that  the  beam  was 
accurately  centered  on  the  detector  at  position  A,  then  rotated  pre¬ 
cisely  90°,  leaving  the  detector  at  position  B. 
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Mirror  Ml  was  then  put  into  place  and  adjusted  until  the  laser 
beam  was  centered  upon  mirror  M2.  Similarly,  M2  was  adjusted  until  the 
beam  was  centered  upon  M3. 

With  the  lenses  LI  and  L2  and  the  pinhole  P  removed  from  the 
system,  mirror  M3  and  mirror  M2  were  systematically  adjusted  until  the 
beam  was  centered  upon  both  the  sample  holder  and  the  detector  at 
position  B. 

The  incident  beam  was  thus  accurately  positioned  perpendicular 
to  the  axis  of  rotation  of  the  movable  arm.  A  polished  sample  was  then 
placed  in  the  holder  and  adjusted  until  the  specularly-reflected  beam 
returned  precisely  along  the  incident  beam.  The  sample  holder  was  thus 
accurately  positioned  perpendicular  to  the  incident  beam  and  locked 
into  this  position. 

Lenses  LI  and  L2  were  then  placed  in  the  beam  and  properly  cen¬ 
tered.  And  finally,  the  pinhole  P  was  accurately  positioned  at  the 
back  focal  position  of  lens  LI. 

With  the  system  properly  aligned,  the  movable  arm  could  be 
rotated  to  direct  the  incident  beam  at  any  desired  angle  without 
requiring  ucher  adjustments  to  keep  the  beam  centered  upon  the  sample. 

An  additional  requirement  was  that  the  PMT  with  its  associated 
fiber-optic  probe  be  positioned  and  aligned  such  that  the  field  of  view 
of  the  PMT  remain  accurately  centered  on  the  illuminated  portion  of  the 
sample  throughout  the  entire  range  of  its  scanning  motion.  Provisions 
were  therefore  made  to  allow  three  degrees  of  freedom  (one  translation 
and  two  rotation)  in  adjusting  the  position  and  orientation  of  the  PMT 
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housing.  Removing  the  PMT  from  its  housing  and  illuminating  the  fiber¬ 
optic  probe  from  the  back  side  greatly  facilitated  this  alignment  pro¬ 
cedure  as  it  allowed  one  to  directly  observe  the  field  of  view  on  the 
sample  holder  while  making  the  necessary  adjustments. 

Measurement  Technique 

The  sample  to  be  measured  was  placed  in  the  holder  and  the  mov¬ 
able  arm  positioned  to  achieve  the  desired  angle  of  incidence.  A  cali¬ 
brated  attenuator  was  then  placed  in  the  incident  beam  and  the  detector 
centered  on  the  specularly-reflected  beam.  The  collecting  aperture  of 
the  fiber-optic  probe  was  large  enough  to  collect  the  entire  specular 
beam;  hence,  the  output  signal,  V8,  of  the  PMT  in  this  position  was 
proportional  to  the  total  flux  in  the  specular  beam. 

The  detector  was  then  moved  a  known  angular  distance  (approxi¬ 
mately  one  degree)  from  the  specular  beam  and  the  attenuator  removed. 

A  profile  of  the  scattered  light  distribution  was  then  measured  by 
scanning  the  observation  hemisphere  with  the  fiber-optic  probe. 

Approximately  30  separate  readings  were  taken  at  different  angular 
% 

positions  between  the  specular  beam  and  the  plane  of  the  sample.  These 
readings  constitute  the  raw  data. 

The  sample  was  then  removed  and  the  incident  beam  allowed  to 
pass  unobstructed  through  the  sample  holder  and  into  a  black  absorbing 
Rayleigh  horn.  Background  measurements  were  then  made  along  the  s.  ^e 
profile  as  above  and  subtracted  from  the  raw  data.  These  background 
measurements  were  found  to  be  completely  negligible  in  most  instances. 
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The  data  now  represents  the  spread  function  of  the  scattering 
system,  which  is  made  up  of  the  spread  function  of  the  scattering 
surface  convolved  with  the  spread  function  of  the  optical  system  pro¬ 
ducing  the  incident  beam.  These  are  shown  in  Fig.  4.11. 

The  spread  function  of  the  incident  beam  is  then  measured  by 
again  placing  the  calibrated  attenuator  in  the  incident  beam  and  cen¬ 
tering  the  detector  on  the  direct  beam  passing  through  the  empty  sample 
holder.  Since  the  collecting  aperture  of  the  fiber-optic  probe  is 

large  enough  to  collect  the  entire  incident  beam,  the  output  signal, 

/ 

V  ,  of  the  PMT  in  this  position  is  proportional  to  the  total  flux  in 
the  incident  beam.  The  detector  is  then  moved  a  known  angular  dis¬ 
tance  (approximately  one  degree)  from  the  direct  beam  and  the  atten¬ 
uator  removed.  A  profile  of  tne  incident  beam  is  then  measured.  These 
readings  rapidly  diminish  to  zero  within  five  degrees  of  the  peak  value. 

Since  the  spread  function  of  the  incident  beam  [Fig.  4.11(a)]  is 
narrow  compared  to  the  scattering  function  of  the  surface  [Fig.  4.11(b)], 
the  scattered  portion  of  the  surface  spread  function  is  virtually  un¬ 
changed  by  the  convolution  operation  while  the  delta  function  component 
merely  replicates  the  beam  spread  function.  The  desired  scattering 
function  can  thus  be  obtained  by  subtracting  the  beam  spread  function 
readings  from  the  raw  data. 

Presentation  of  the  Data 

It  is  customary  to  present  scattered  light  data  from  diffusely 
reflecting  surfaces  in  a  polar  format.  Three  experimental  curves  and 
one  ideal  reference  curve  are  illustrated  in  Fig.  4.12.  Note  that  we 
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Fig.  4.11.  Illustration  of  Components  Comprising 
System  Spread  Function. 

(a)  Spread  Function  of  Incident  Beam. 

(b)  Spread  Function  of  Scattering  Surface. 

(c)  Spread  Function  of  Scattering  System. 


have  plotted  the  quantity  y — —  along  the  radial  coordinates  whore 

o 

Au'  is  the  solid  angle  subtended  by  the  fiber-optic  probe.  Since  V  is 
proportional  to  the  power  collected  by  the  fiber-optic  probe  and  V 
is  proportional  to  the  total  power  in  the  incident  beam,  this  quantity 
is  the  relative  intensity  cf  the  scattered  light  distribution  (i.e., 
scattered  intensity  normalized  by  the  incident  power).  For  smooth 
mirror  surfaces  this  same  quantity  is  usually  plotted  in  a  cartesian 
format  as  a  function  of  the  scattering  angle  as  shown  in  Fig.  4.13. 

Dividing  this  data  by  the  cosine  of  the  scattering  angle 
(y  a  cos0) ,  we  obtain 


V  P/  G4yAui) 

V0L  “Y  ~  Pq/A  ~  Ei 


(73) 


where  A  is  the  illuminated  area  on  the  sample.  This  quantity  is 
equivalent  to  the  reflected  radiance  in  the  sampled  direction  divided 
by  the  incident  irradiance,  which  is  precisely  the  manner  in  which  the 
BRDF  is  defined.  The  resulting  scattering  curve,  exhibited  in  Fig. 

4.14,  is  therefore  a  one-dimensional  profile  of  the  four-dimensional 
BRDF.  Also,  in  accordance  with  our  theory  we  are  plotting  this  function 
versus  the  quantity  8-8  >  which  is  the  distance  of  the  observation 
point  from  ihe  specular  beam  in  direction  cosine  space. 

It  can  be  readily  shown  that  the  BRDF  is  merely  an  infinite 
family  of  two-dimensional  spread  functions  which  are  scaled  by  the 
total  reflectance  of  the  surface. 
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Fig.  4.13.  Illustration  of  Relative  Intensity  Plotted 
as  a  Function  of  Scattering  Angle. 
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Fig.  4.14.  Illustration  of  BRDF  Profiles 
Plotted  as  a  Function  of  B-3o- 
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BRDF 


(74) 


=  X  7?-£(ct,B), 


Since  we  measure  only  the  scattered  component  of  the  spread  function 
(See  Eq.  57)  and  are  primarily  concerned  with  the  scattered  light 
behavior  rather  than  the  total  reflectance,  the  scattering  function 
can  be  written  as 


S(ct,6)  = 


R  V0Auj  y 


All  measurements  are  thus  normalized  by  the  reflectance  of  the  surface 
so  it  does  not  appear  to  have  better  scattering  characteristics  due  to 
its  lower  reflectance. 

The  MgO  surface,  which  is  a  fairly  good  diffuse  reflectance 
standard,  yields  a  straight  horizontal  line  as  a  Lambertian  reflector 
should.  Since  the  radiance  of  a  Lambertian  source  is  given  by 


where  M  is  the  total  emittance  into  a  hemisphere,  the  value  of  the 
MgO  scattering  function  can  be  shown  to  be  ^  .  A  Lambertian  surface 
of  known  reflectance,  R £,  thus  makes  a  convenient  reference  sample 
and  the  scattering  function  of  a  test  sample  of  reflectance,  R,  is 
given  by 

,  7(a,8)  RL 

S(  a,6)  *  - L  (77 

R  VL  n  y 

where  V £  is  trie  PMT  voltage  signal  from  the  Lambertian  reference  sur¬ 
face  at  a=g=0.  A  freshly  coated  MgO  sample  with  a  total 
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hemispherical  reflectance  J?£  =  0.98  was  routinely  used  as  a  reference 
sample  for  the  remainder  of  the  measurements.  Unless  stated  otherwise, 
the  scattering  function  of  Eq.  (77)  will  be  used  for  presenting  the 
data  in  this  dissertation.  Furthermore,  in  accordance  with  the  theory 
presented  in  Chapter  2,  this  scattering  function  will  be  plotted  as  a 
function  of  distance  from  the  specular  beam  in  direction  cosine  space 
(8-B0)  as  shown  in  Fig.  4.15.  Note  that  by  plotting  both  ordinate  and 
abscissa  on  a  log  scale  the  scattering  function  for  the  polished  sample 
takes  the  form  of  a  straight  line. 
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Scattering  Function 


CHAPTER  5 


DATA  AND  RESULTS 

The  angular  distribution  of  scattered  light  measured  from  a 
variety  of  surfaces,  angles  of  incidence,  and  wavelengths  are  reported 
in  this  chapter.  The  results  obtained  are  then  compared  with  theoreti¬ 
cal  predictions. 


Surface  Roughness  Effects 

The  scattered  light  profile  for  a  normally  incident  beam  on  a 
variety  of  samples  with  a  wide  range  of  rms  surface  roughness  values 
are  exhibited  in  Fig.  5.1.  Note  that  the  polished  samples  are  charac¬ 
terized  by  a  straight  line  curve  with  a  slope  between  -3/2  and  -2.  The 
ground  glass  samples  yield  curves  which  are  flat  for  a  substantial 
angular  range  before  falling  off  rapidly  at  the  larger  angles.  We 
were  unable  to  produce  samples  that  satisfactorily  bridged  the  obvious 
gap  in  the  data  between  the  ground  and  polished  samples.  A  separate 
study  of  the  grinding  and  polishing  process  on  fused  silica  (Shevlin, 
1974)  utilizing  electron  microscopic  examination  of  the  surface  also 
indicates  a  very  rapid  change  in  the  surface  character  between  the  fine 
grind  and  the  polishing  operation.  This  is  dramatically  illustrated  in 


xerographic  Display  Illustrating  the  Surface  Structure  of  Fused 
Various  Stages  of  the  Grinding  and  Polishing  Process. 


It  is  also  instructive  to  compare  the  scattering  curves 
resulting  from  special  materials  or  unusual  fabrication  techniques  with 
those  of  more  conventional  optical  surfaces.  For  example,  the  results 
of  scattering  measurements  on  a  polished  beryllium  sample  and  an 
ionically-polished  fused  quartz  sample  are  shown  in  Fig.  5.3  along  with 
some  curves  from  conventional  optical  surfaces. 

Incident  Angle  Effects 

Rather  extensive  scatter  measurements  have  been  made  on  two 
representative  surfaces.  One  is  a  ground  glass  surface  (Sample  #172) 
which  is  a  very  diffuse  reflector.  The  other  is  an  optically-polished 
surface  (Sample  #200)  which  is  a  nice  specular  reflector.  Both  samples 
were  coated  with  aluminum  prior  to  making  the  measurements.  Four 
separate  scattering  profiles  from  the  specular  beam  to  the  plane  of  the 
sample  (see  Fig.  4.7)  were  measured  at  several  angles  of  incidence. 

The  backscattering  profile  of  the  scattered  light  distribution 
for  these  two  samples  is  shown  in  Fig.  5.4  for  several  different  angles 
of  incidence.  For  the  polished  sample,  the  various  curves  coincide 
almost  perfectly  for  angles  of  incidence  between  zero  and  60°.  Hence, 
it  is  apparent  that  the  scattering  function  does  not  appreciably  change 
with  the  angle  of  incidence.  The  corresponding  curves  for  the  rough 
sample  coincide  for  a  substantial  range  of  angles  then  begin  to  depart 
somewhat  at  the  large  angles. 

The  four  separate  profiles  of  the  scattered  light  distribution 
fro..  the  same  two  samples  with  the  incident  beam  at  45°  are  shown  in 
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Scattering  Function 


•  Conventional  Polish 

A  Ion  polished 
C  Polished  beryllium 
^  Submerged  polish 
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Tig.  5.3.  Comparison  of  Scatter  Profiles  from  Special  Materials 
and  Unusual  Fabrication  Techniques  with  Those  of 
Conventional  Optical  Surfaces. 
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Scattering  Function 


Fig.  5.4.  Comparison  of  Scatter  Profiles  for  Different  Incident  Angl 
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Fig.  5.5.  Again  the  curves  for  the  polished  sample  coincide  almost 
perfectly,  suggesting  a  rotationally-symmetric  distribution  in  direc¬ 
tion  cosine  space.  Some  asymmetry  is  noted  in  the  scattered  light 
distribution  from  the  rough  sample. 

The  data  on  Figs.  5.4  and  5.5  confirm  that  for  a  certain  class 
of  surfaces  (in  which  optically-polished  glass  is  definitely  a  member, 
and  ground  glass  can  perhaps  be  included  to  a  lesser  extent),  the  scat¬ 
tering  properties  are  indeed  shift  invariant  as  predicted  by  the  theory 
presented  in  Chapter  3,  and  can  be  completely  characterized  by  a  single 
set  of  measurements  at  a  fixed  angle  of  incidence! 

Figure  5.6  dramatically  illustrates  the  importance  of  the 
coordinate  system  within  which  the  scattering  process  is  discussed. 

The  curves  in  Fig.  5.6(a)  correspond  to  the  scattered  light  distribution 
illustrated  in  Fig.  3.4  and  confirms  the  well-known  fact  that  a  curve 
representing  the  scattered  intensity  as  a  function  of  observation  angle 
will  change  shape  drastically  with  angle  of  incidence--becoming  quite 
skewed  and  asymmetrical  at  large  angles  of  incidence.  However,  these 
same  data,  when  plotted  in  accordance  with  our  theory,  describes  a  new 
scattering  function  which  does  not  change  shape  but  will  merely  be 
shifted  in  direction  cosine  space  with  changes  in  angle  of  incidence  as 
shown  in  Fig.  5.6(b).  This  is  a  rather  significant  development  which 
greatly  reduces  the  quantity  of  data  required  to  completely  character¬ 
ize  the  scattering  properties  of  a  surface. 
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Fig.  5.5.  Comparison  of  Scatter  Profiles  Taken  in  Different 
Directions  for  an  Incident  Beam  at  45  Degrees. 


The  Inverse  Scattering  Problem 


The  problem  of  determining  surface  characteristics  from  scat¬ 
tered  light  measurements  plays  an  important  role  in  many  areas  of 
technical  interest.  Recall  that  the  surface  autocovariance  function 
and  the  rms  surface  roughness  are  two  relevant  surface  characteristics 
which  are  related  to  the  scattering  properties  of  the  surface  by  the 
transfer  function  described  in  Chapter  3. 

Predictions  of  Surface  Characteristics 

A  computer  program  has  been  written  (see  Appendix  C)  for  cal¬ 
culating  the  effective  transfer  function  and  the  surface  autocovariance 
function  from  scattered  light  data.  The  measured  data  are  assumed  to 
be  a  radial  profile,  S(p),  of  a  rotationally-symmetric  scattering  func¬ 
tion.  An  intermediate  quantity,  BQ(8),  is  first  determined  by  calcula¬ 
ting  the  two-dimensional  Fourier-Bessel  Transform  of  this  scattering 
function.  The  surface  transfer  function  is  then  calculated  from 
Eq.  (53),  where  A  =  1 -B.  And  finally,  the  surface  autocovariance  func¬ 
tion  is  computed  from  Eq.  (58).  Representative  curves  for  a  diffusely 
reflecting  ground  glass  surface  and  a  specularly  reflecting  polished 
glass  surface  are  shown  in  Fig.  5.7  and  Fig.  5.8.  Note  that  the  total 
integrated  scatter,  B,  and  the  rms  surface  roughness,  o^,  are  also  pro¬ 
vided  from  these  curves.  Figures  5.9  through  5.13  illustrate  the  pre¬ 
dicted  surface  characteristics  from  other  samples  which  we  wish  to 
compare  with  independent  measurements. 
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Surface  Autocovariance 
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Surface  Autocovariance  Function  (*  10-6  yin 


Surface  Autocovariance  Function 


Surface  Structure  Measurements 


Three  traditional  techniques  for  obtaining  surface  structure 
information  involve  profilometry,  electro-microscopy,  and  FECO  inter¬ 
ferometry.  Most  profi ’ometers  provide  too  coarse  a  measurement  for 
optical  surfaces.  As  illustrated  in  Fig.  5.2  the  electron  microscope 
works  nicely  on  the  rough  ground  glass  surfaces  but  fails  to  yield 
sufficient  information  about  the  smooth  polished  surfaces.  The  FECO 
interferometer  works  well  on  smooth  surfaces  with  a  strong  specular 
beam  but  does  not  yield  good  results  for  the  rough  diffusely  reflecting 
surfaces  since  it  requires  multiple  reflections.  The  latter  two  com¬ 
plementary  techniques  were  thus  utilized  in  our  research  effort. 

Surface  profiles  of  rough  samples  can  be  determined  from 
electron-micrograph  stereo  pairs  using  conventional  stereo- 
photogrammetric  techniques  (Moffitt,  1959).  Nankivell  (1963)  discusses 
some  of  the  stereo-photogrammetric  problems  unique  to  electron  micro¬ 
scope  applications.  Electron-micrograph  stereo  pairs  were  produced 
with  both  a  conventional  transmission  electron  microscope  and  a 
scanning  electron  microscope  at  a  variety  of  tilt  angles  and  magnifica¬ 
tions  for  several  surfaces  with  known  characteristics.  Considerable 
preliminary  experimentation  with  stereo-photogrammetric  techniques  was 
then  performed  with  a  variety  of  stereoscope-parallax  bar  configurations. 

A  typical  electron-micrograph  stereo  pair  is  shown  in  Fig.  5.14 
with  a  line  scribed  to  indicate  the  position  of  a  set  of  preliminary 
surface  height  measurements  that  were  made  with  a  standard  Fairchild 
Stereocomparagraph.  This  instrument  consists  of  a  mirror  stereoscope 
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fitted  with  a  parallax  bar  containing  a  micrometer  for  measuring  the 
parallax  of  each  desired  pair  of  points.  The  resulting  surface  profile 
is  also  shown  in  Fig.  5.14.  This  surface  is  a  very  rough  ground  glass 
produced  with  30  nm  grit.  The  tedious  procedure  of  obtaining  surface 
profiles  in  this  manner  becomes  increasingly  more  difficult  as  the 
surface  becomes  smoother. 

The  two  statistical  parameters  which  determine  the  scattered 
light  characteristics  are  the  variance  of  the  surface  height  distribu¬ 
tion  and  the  surface  autocorrelation  function.  A  computer  program  (see 
Appendix  D)  was  written  which  takes  the  surface  profile  data  and 
determines  the  above  two  parameters.  An  electron  micrograph  showing 
the  surface  of  Sample  #172  magnified  10,000  times  is  illustrated  in 
Fig.  5.15  along  with  the  surface  profile,  surface  height  distribution 
function,  and  the  surface  autocovariance  function.  One  hundred -twenty 
data  points  were  used  for  obtaining  these  surface  parameters.  Addi¬ 
tional  data  was  recorded  to  determine  the  degree  to  which  the  surface 
was  homogeneous  and  isotropic. 

Surface  profiles  of  several  smooth  samples  were  measured  by 
Dr.  Jean  Bennett  at  the  Naval  Weapons  Center  on  a  scanning  FECO  inter¬ 
ferometer  (see  Appendix  E)  capable  of  determining  very  small  height 
differences  with  a  lateral  resolution  of  2  ym.  This  instrument,  along 
with  auxiliary  equipment  which  includes  a  slow-scan  TV  camera,  signal 
averager,  minicomputer,  and  teletype  unit,  yields  the  surface  profile, 
rms  roughness,  surface  height  distribution  function,  surface  autocovar¬ 
iance  function,  and  other  statistical  properties  of  optical  surfaces. 
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Electron  Micrograph 


An  Electron  Micrograph  of  Sample  #172  is  shown  along  with  a  Surface 
Profile,  the  Surface  Height  Distribution,  and  the  Surface  Autocovar 
ance  Function  as  Determined  by  Stereo-Photogrammetric  Techniques. 
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The  results  of  these  measurements  on  Sample  #198  are  shown  in 
Fig.  5.16. 

The  rms  surface  roughness  was  also  predicted  from  visual  mea¬ 
surements  of  the  FECO  fringe  widths,  as  well  as  from  total  integrated 
scatter  measurements  performed  at  the  Naval  Weapons  Center.  Additional 
scattered  light  measurements  made  on  a  few  samples  allowed  the  autoco¬ 
variance  length  to  be  calculated  when  assuming  the  autocovariance 
function  to  be  Gaussian. 


Comparison  of  Predictions  with  Measurements 

Theoretical  predictions  of  surface  characteristics  for  Sample 
#172  were  shown  in  Fig.  5.7.  Direct  comparison  with  the  results  of  the 
surface  measurements  shown  in  Fig.  5.15  indicate  that  the  values  for 
the  rms  surface  roughness  differ  by  approximately  a  factor  of  three  and 
the  widths  of  the  autocovariance  function  (distance  to  the  first  zero) 
differ  by  approximately  20  per  cent.  Under  the  circumstances  involved 
this  can  be  considered  to  be  remarkably  good  agreement  since  this  sur¬ 
face  is  so  rough  as  to  make  the  theory  somewhat  suspect  and  not  rough 
enough  for  the  electron-micrograph  stereo  measurements  to  be  considered 
reliable. 

The  predictions  of  surface  characteristics  based  upon  scattered 
light  measurements  from  smooth  surfaces  are  tabulated  in  Table  5.1  along 
with  the  results  of  independent  measurements  performed  at  the  Naval 
Weapons  Center.  The  rms  surface  roughness  of  a  few  samples  was 
determined  both  from  total  integrated  scatter  measurements  and  from 
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Fig.  5.16, 


Surface  of  Sample  #198  as  Determined  by  the  Scanning 
FECO  Interferometer  at  the  Naval  Weapons  Center. 
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Table  5.1.  Surface  Characteristics. 
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visual  measurement  of  the  PECO  fringe  widths.  The  wide  variation  in 
the  results  of  these  two  techniques  for  Sample  #186  and  Sample  #163  was 
attributed  to  nonuniform  roughness  over  the  surface  of  the  sample.  The 
surface  autocovariance  width  determined  from  the  scanning  FECO  inter¬ 
ferometer  is  one  to  two  orders  of  magnitude  larger  than  that  theoreti¬ 
cally  predicted  from  scattered  light  measurements.  This  is  due  to 
the  poor  lateral  resolution  of  the  instrument.  The  surface  autocovari¬ 
ance  width  predicted  from  a  simple  theory  which  assumes  a  Gaussian 
shape  for  the  autocovariance  function  does  compare  favorably  with  our 
predictions.  The  autocovariance  widths  tabulated  in  Table  5.1  are 
arbitrarily  chosen  to  be  the  half-width  of  the  surface  autucovariance 
function  at  l//e  times  its  maximum  height. 

Indirect  Verification  of  Inverse  Scattering  Solution 

The  accurate  determination  of  micro- structure  on  an  optical 
surface  has  been  shown  to  be  extremely  difficult  to  accomplish  by 
direct  measurement  (therein  lies  the  motivation  for  attempting  to  solve 
the  inverse  scattering  problem) .  However,  our  theoretical  treatment  of 
the  inverse  scattering  problem  can  be  tested  indirectly  by  using  the 
surface  characteristics  predicted  from  scattered  light  measurements  at 
a  given  wavelength  to  calculate  the  scattered  light  behavior  at  a  dif¬ 
ferent  wavelength.  This  calculated  scattering  function  can  then  be 
readily  compared  to  the  directly  measured  scattering  function  at  that 
wavelength.  Figure  5.17  provides  a  direct  comparison  of  the  measured 
scattering  function  from  Sample  #184  at  X  =  0.6328  ym  to  that  predicted 
from  surface  characteristics  determined  from  scattered  light 
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Sample  #184 


measurements  at  X  =  0.4579  ym.  The  remarkable  agreement  displayed  pro¬ 
vides  experimental  verification  of  our  theoretical  treatment  of  the 
inverse  scattering  problem,  at  least  for  smooth  surfaces. 

Wavelength  Effects 

The  transfer  function  characterization  of  scattering  surfaces 
developed  in  Chapter  3  offers  a  simple  means  of  determining  the  wave¬ 
length  dependence  of  the  scattered  light  behavior.  Once  the  surface 
characteristics  are  known  (whether  from  direct  measurement  or  calcu¬ 
lated  from  scattered  light  data),  the  same  theory  provides  a  simple 
method  of  predicting  the  scattering  function  at  any  desired  wavelength. 

This  technique  of  predicting  the  scattering  behavior  as  a  func¬ 
tion  of  wavelength  has  been  experimentally  verified  (see  Fig.  5.17)  for 
wavelength  ranges  limited  to  the  visible  spectrum.  Furthermore,  the 
same  data  can  be  used  to  verify  the  wavelength  scaling  law  for  smooth 
surfaces  stated  in  Eq.  (70).  This  is  illustrated  in  Fig.  5.18. 

Similar  attempts  to  predict  the  scattering  properties  of  a  sur¬ 
face  at  a  wavelength  of  10.6  ym  from  measured  data  in  the  visible  was 
not  successful.  This  failure  was  due  to  the  greatly  expanded  angular 
width  of  the  scattering  function  at  long  wavelengths  as  described  by 
the  wavelength  scaling  law,  and  is  illustrated  in  Fig.  5.19.  Note  that 
scattered  light  measurements  over  angular  range  from  1°  to  46°  at  a 
wavelength  of  0.4579  ym  can  be  used  to  predict  the  scattering  behavior 
from  1.4°  to  90°  for  a  wavelength  of  0.6328  ym.  However,  these  same 
measurements  provide  only  a  few  data  points  in  the  angular  range  from 
1°  to  2.5°  that  are  useful  in  predicting  the  scattering  behavior  for  a 


111 


Sample  #184 
6n  =  0 


wavelength  of  10.6  pm,  and  no  information  is  obtained  concerning  the 
scattered  light  behavior  at  angles  less  than  24°. 

It  is  now  clear  that  it  may  not  be  possible  to  compare  the 
scattered  intensity  of  two  widely  separated  wavelengths  in  a  given 
direction  without  extiapolating  one  curve.  However,  this  behavior  has 
the  advantage  of  allowing  one  to  determine  the  scattering  characteris¬ 
tics  at  very  small  angles  (unobtainable  by  direct  measurement  due  to 
mechanical  constraints)  by  making  large  angle  scatter  measurements  at  a 
longer  wavelength.  The  angular  range  of  validity  for  predictions  of 
scattered  light  behavior  based  upon  measurements  from  1°  to  80°  at 
10.6  pm  is  presented  as  a  function  of  wavelength  in  Fig.  5.20.  For 
example,  measurements  from  1°  to  80"  at  a  wavelength  of  10.6  pm  could 
be  used  to  predict  the  scattering  properties  of  visible  light 
(A  =  0.5  pm)  in  an  angular  range  from  approximately  0.04o^  to  2.7°.  Or 
conversely,  if  one  has  the  capability  of  making  very  small  angle  mea¬ 
surements  in  the  visible,  the  wide  angle  scattering  characteristics  at 
a  longer  wavelength  can  be  determined. 


Scattering  Angle 


Fig.  5.20.  The  Angular  Range  of  Validity  for  Predictions  of 
Scattered  Light  Behavior  Based  upon  Measurements 
from  1°  to  80°  at  10.6  ym. 
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CHAPTER  6 

CONCLUSION 

This  study  has  been  a  general  investigation  of  surface  scatter 
phenomena  dealing  with  several  different  aspects  of  scattered  light 
behavior.  An  elementary  theoretical  development  based  upon  scalar  dif¬ 
fraction  theory  has  been  presented.  Linear  systems  theory  and  modern 
Fourier  techniques  result  in  a  theoretical  model  of  light  scattering 
systems  which  closely  parallels  the  highly  successful  theory  of  iso- 
planatic  imaging  systems.  An  extensive  experimental  program  has 
accompanied  this  theoretical  development  in  an  attempt  to  verify 
theoret i ca 1  predict ions. 

Summary  of  Results 

An  analytical  expression  has  been  obtained  for  a  surface  trans¬ 
fer  function  which  relates  the  surface  micro-structure  to  the  scattered 
distribution  of  radiation  from  that  surface.  The  existence  of  such  a 
transfer  function  implies  a  shift-invariant  scattering  function  which 
does  not  change  shape  with  the  angle  of  the  incident  beam.  This  is  a 
rather  significant  development  which  greatly  reduces  the  quantity  of 
data  required  to  completely  characterize  the  scattering  properties  of  a 
surface.  For  a  large  class  of  well-behaved  surfaces  this  transfer 
function  is  described  in  term.'  of  only  the  rms  surface  roughness  and 
the  surface  autocovariance  function.  This  transfer  function  thus 
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provides  a  straightforward  solution  to  the  inverse  scattering  problem 
(i.e.,  determining  the  surface  characteristics  from  scattered  light 
measurements).  Once  the  surface  characteristics  are  known,  the  same 
theory  provides  an  equally  simple  method  of  predicting  the  wavelength 
dependence  of  the  scattered  light  distribution. 

Experimental  verification  of  the  shift-invariant  scattering 
function  has  been  successfully  demonstrated  for  smooth  surfaces  (o^«X) 
The  scattered  light  measurements  from  rough  (diffusely  reflecting)  sur¬ 
faces  results  in  a  scattering  function  which  is  shift- invariant  over 
only  a  small  range  of  angles  and  departs  significantly  from  the  pre¬ 
dicted  behavior  at  large  scattering  angles. 

A  computer  program  has  been  developed  that  operates  upon 
scattered  light  data  to  yield  the  total  integrated  scatter,  the  surface 
transfer  function,  the  rms  surface  roughness,  and  the  surface  autoco¬ 
variance  function.  Although  accurate  determination  of  micro-structure 
on  optical  surfaces  is  extremely  difficult  to  accomplish  by  direct 
measurement  (thus  the  motivation  for  attempting  to  solve  the  inverse 
scattering  problem),  favorable  comparisons  of  predicted  surface  charac¬ 
teristics  with  the  corresponding  measured  quantities  have  been  demon¬ 
strated  for  both  smooth  surfaces  and  moderately  rough  surfaces.  In 
addition,  experimental  verification  of  the  inverse  scattering  program 
was  accomplished  indirectly  by  supplying  scattered  light  data  of  one 
wavelength  as  input  to  the  inverse  scattering  program  in  order  to  deter 
mine  the  relevant  surface  characteristics;  then  this  information  was 
used  to  predict  the  scattering  function  at  a  different  wavelength. 
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Excellent  agreement  with  the  measured  scattering  function  at  that 
wavelength  was  achieved. 

Since  the  above  technique  involves  numerical  computations  on 
sampled  data,  an  analytical  expression  for  a  wavelength  scaling  law  is 
not  required  to  determine  the  scattering  function  at  any  desired  wave¬ 
length.  However,  in  order  to  gain  insight  into  the  wavelength  depen¬ 
dence  of  surface  scatter  phenomena,  a  wavelength  scaling  law  for  smooth 
surfaces  was  derived  and  verified.  This  scaling  law  consists  of  a 
change  in  the  scattering  angle  as  well  as  a  change  in  the  amplitude  of 
the  scattering  function  with  changes  in  wavelength.  It  therefore  pro¬ 
vides  a  valuable  tool  for  predicting  the  scattering  behavior  in  certain 
angular  regions  or  wavelength  ranges  where  direct  measurements  are 
difficult  to  obtain. 

Further  Theoretical  Considerations 

Ho  explicit  approximations  concerning  the  size  of  the  surface 
variations  were  made  in  the  theoretical  development  presented  in  Chap¬ 
ter  3.  However,  a  simplifying  assumption  was  made  regarding  the  random 
component  of  the  pupil  function  described  in  Eq.  (33).  It  was  assumed 
that  the  phase  variations  in  the  disturbance  emerging  from  the  scatter¬ 
ing  surface  were  equal  to  the  perturbations  introduced  onto  a  normally 
incident  wavefront.  Careful  examination  of  Fig.  6.1  reveals  that  the 
phase  difference  introduced  by  reflection  from  a  rough  surface  depends 
upon  both  the  angle  of  incidence  and  the  angle  of  observation  in  addi¬ 
tion  to  the  surface  height  at  the  point  of  reflection.  The  phase 
variations  along  the  scattering  surface  can  thus  be  expressed  as 


t/ 


Surface 


Fig.  6.1.  Illustration  of  the  Phase  Variation 
Introduced  by  Reflection  from  a 
Rough  Surface. 
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where 


4 >(x,y)  =  2tt(y+yo)Wb(£,&) 


Y  =  cose,  y  Q  =  coseQ. 


(78) 


The  effective  transfer  function  of  a  scattering  surface  is  now 
given  by  the  following  general  expression 


'X’xQp.y)  c  e 


-[2*(y+v>av  2] 


!  .rnMl 


y 


C79) 


This  expression  can  be  interpreted  as  a  two-parameter  family  of  trans¬ 
fer  functions,  one  for  every  possible  angle  of  incidence  and  every 
possible  scattering  angle. 

This  generalization  still  leaves  us  with  a  theoretical  model 
closely  paralleling  that  of  non-isoplanatic  imaging  systems  which  can 
be  characterized  by  a  different  transfer  function  for  each  off- 
axis  object  point.  It  therefore  seems  reasonable  that  an  "aberration 
theory"  of  scattering  systems  can  be  developed  to  provide  more  insight 
into  the  scattering  behavior  of  rough  surfaces. 

A  preliminary  empirical  search  for  a  new  scattering  function 
of  the  form  (y+yq)  ,8)  plotted  versus  the  quantity  (8-8o)/(Y+Y0)n 
resulted  in  the  following  interesting  dilemma.  When  m= 2  and  n= 1,  the 
•  scattering  function  for  Sample  #172  indeed  becomes  quite  shift-invariant 
with  respect  to  changes  in  the  incident  angle  (there  is  a  theoretical 
basis  for  this  particular  form  of  the  scattering  function  if  some 
assumptions  are  made  concerning  the  shape  of  the  surface  autocovariance 
function) .  However,  this  same  scattering  function  degrades  the 


shift-invariance  of  Sample  #200.  Clearly  a  more  general  theory  should 
also  work  for  the  special  case  of  smooth  surfaces. 

Suggestions  for  Further  Research 
The  following  suggestions  are  made  for  further  research  in  the 
area  of  surface  scatter  phenomena: 

1.  Continue  the  theoretical  development  on  the  transfer  function 
characterization  of  scattering  surfaces.  This  should  include  an 
attempt  to  solve  the  inverse  scattering  problem  for  rough  surfaces 

by  using  the  general  expression  for  the  transfer  function  described 
by  Eq.  (79).  Empirical  curve  fitting  techniques  may  be  useful  in  an 
attempt  to  discover  new  scattering  functions  or  plotting  techniques 
which  result  in  shift-invariant  behavior  with  respect  to  changes  in 
the  incident  angle.  Generalizing  the  surface  scatter  theory  to  a 
complete  vector  treatment  while  maintaining  the  transfer  function 
approach  would  be  a  major  contribution  to  the  understanding  of  sur¬ 
face  scatter  phenomena. 

2.  Improve  the  scattered  light  measuring  capabilities  by 
obtaining  laser  sources  that  span  a  larger  wavelength  range. 

3.  Automate  the  scattering  apparatus  for  high-speed  acquisi¬ 
tion  and  analysis. 

4.  Acquire  samples  and  perform  measurements  upon  selected 
moderately  rough  surfaces  (0.1<a<1.0)  and  compare  with  theoretical 
predictions. 


121 


APPENDIX  A 


CALCULATION  OF  ABERRATION  COEFFICIENTS 
FOR  DIFFRACTED  WAVE  FIELDS 


For  the  case  of  a  plane  wave  incident  upon  an  aperture,  the 
diffracted  wave  field  on  an  observation  plane  is  given  by  Eqs.  (16)  and 
(17),  where 

W  «  (£-&)  +  (£'£  +  $'£)/2.  (A. 1) 

* 

The  quantity  l  can  be  written  as 

t  =  >/(£-£')  2  +  (p-p')2  + 


=  6  A  +  -  2I5F+PPTI7I2',  (A.  2) 


where 


§2  =  &2  +  g2 ,  gi2  -  £t2  +  p'2t 

A  binomial  expansion  of  the  above  square  root  results  in  the  following 
expression  for  W 

W  =  |  [§2  +  §'2]/S2 

-  |  [$k+$'l*+4{&£'+gg')2+2$2$'2-te2{£&'+gQ')-4s'2(xx'+yy')]/2k 
+  higher-order  terms.  (A. 3) 
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If  we  assume  a  rotational ly- symmetric  diffracting  aperture  we 
can,  without  loss  of  generality,  choose  the  observation  point  on  the 
p-axis.  Let  us  therefore  set  x  =  0.  We  can  then  let  s  =  p  and 
p’  a  S'  cos*,  hence 


W 


|  [y2  ♦  s'2]/a2 


Jr  IP4  +  s'4  +  4p2a'2  cos2<{>  +  2p28,2-4p3s'  cos4> 


+  higher-order  terms. 


-  4p§'3  cos*] /a4 

O'  *4) 


If  we  now  substitute 


A  A  A  |  A  Ct 

V  *  p^max’  8  B 


intc  the  previous  equation,  we  obtain 


i/  =  |  [sBax2e2*  e/2)2  a2]/®2 


-  |  [S„ax4  P-  *  (3/2)“  a"  *  45^2  (5/2) 2  o2a2  cos2* 


'max 


+  2pmax2  (3/2)2  p2a2  -  4Pmax3 (5/^)  P3a  cos* 


’  p^3  cos$]/s4 

lUclA 


+  higher-order  terms. 


(A.  5) 


Equating  coefficients  of  corresponding  terms  between  this  equa¬ 
tion  and  the  wavefront  aberration  function  given  by  Eq.  (24),  we  obtain 
the  aberration  coefficients  tabulated  in  the  first  column  of  Table  2.1. 
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If  we  now  have  a  spherical  wave  incident  upon  the  aperture,  the 


quantity  W  in  Eq.  (17)  is  given  by 

W  =  (£-2)  -  (t0-z)  +  (&&'  +  yy')/z. 


(A.  6) 


where 


/£' 2  +  p» 2  +  §2  =  §  /I  +  g'2/g2. 


A  binomial  expansion  of  this  quantity  results  in 


t0  -  a 


(s'/s)2  -  4  (s'/3)4  +  .. 


(A.  7) 


But  these  will  merely  cancel  identical  terms  in  the  previous  expansion 
thus  resulting  in  the  aberration  coefficients  presented  in  column  two  of 


Table  2.1. 


For  the  case  of  a  hemispherical  observation  space  the  diffracted 
wave  field  is  given  by  Eqs.  (22)  and  (23).  With  a  plane  wave  incident 
upon  the  aperture,  we  have 


w  =  (j e-r)  +  (ax'  +  ep  •) 

The  quantity  l  can  be  written  as 


(A.  8) 


where 


l  =  /(£-£')2  +  'fo -O') 2  +  S2 


r  /l  +  [§'2  -  2 w  +  ($')]/*2  , 


•2  =  x'?-  +  y  '2.  r2  =  x2  +  u2  +  s2 


(A.  9) 


r4-  =  +  y4-  +  s4 


i 


t 
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A  binomial  expansion  of  the  above  square  root  yields 


W  =  4  (s'/*)2  -  7  [s'1*  -  4rs' 2(a£'+8y ')  +  4r2(aic'+Bp')2]/f'4 


+  higher-order  terms. 


(A. 10) 


If  we  again  assume  a  rotational ly-symmetric  diffracting  aperture 
we  can,  without  loss  of  generality,  choose  the  observation  point  on  the 
y-axis.  Let  us  therefore  set  a  *  0.  We  can  also  let  y'  «  s'  cos<f>  which 
results  in 

W  -  j  (s'/*)2  -  ^  [s'4  -  4f*gs'3  cos<j>  +  4r282s'2  cos2<j>]/£4 


+  higher-order  terms. 


(A. 11) 


If  we  now  substitute 


PS 


max' 


s'  =  a 


into  the  previous  equation,  we  obtain 


W 


~  (4/ 2r)2  a2  -  j  [ (d/2)1*  a4  -  4r8max  (d/2)3  6p3  cos<f> 


+  4r2f3  2  (d/2)2  82p2  cos2*]/*4 


+  higher-order  terms. 


(A. 12) 


Again  equating  coefficients  of  corresponding  terms  between  this 
equation  and  the  wavefront  aberration  function  given  by  Eq.  (24),  we 
obtain  the  aberration  coefficients  tabulated  in  column  three  of 
Table  2.1. 
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If  we  now  have  a  spherical  wave  illuminating  the  aperture  and  a 
hemispherical  observation  space,  the  quantity  V  in  Eq.  (23)  is  given  by 

W  =  (l-r)  -  (l0-r)  +  Cobs'  +  #•),  (A.  13) 

where 


A  binomial  expansion  of  this  quantity  results  in 

l0  -  r  ■  |  (S'/r)2  -  |  (e Vf*)4  +  ...  .  (A.14) 

Once  again  these  terms  merely  cancel  identical  terms  in  the  previous 
expansion,  leaving  only  coma  and  astigmatism  present  in  the  diffracted 
wave  field  as  indicated  in  the  last  column  of  Table  2.1. 
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APPENDIX  B 
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BIDIRECTIONAL  REFLECTANCE  DISTRIBUTION  FUNCTION 

The  basic  quantity  that  characterizes  (geometrically)  the 
reflecting  properties  of  a  surface  element  dA  is  the  bidirectional 
reflectance  distribution  function  (BRDF).  This  quantity 

-  ;  0j», (sr  *) 

(B.l) 

is  defined  by  Nicodemus  (1970),  as  the  reflected  radiance 
<&riH>H>  *r*$r»  Eti  oi  the  surface  element  dA  in  the  direction  (0r ,<}>*.) 
divided  by  the  incident  irradiance  dEii^iAi)  *  Hi H»H)d&i  producing 
it.  The  geometry  of  this  situation  is  illustrated  in  Fig.  B.l,  where 
the  element  of  projected  solid  angle  is  given  by  d f>  =  cos 0cfu>. 

The  numerical  value  of  the  BRDF  for  a  given  pair  of  incident  and 
reflected  ray  directions  may  vary  frcm  zero  to  infinity.  In  particular, 
consider  two  ideal  cases.  The  BRDF  is  a  constant  for  all  reflected 
directions  for  a  perfectly  diffuse  (Lambertian)  surface;  and  it  becomes 
infinite  (as  a  Dirac  delta  function)  for  a  perfectly  specular  reflector. 
The  BRDF,  defined  above  as  a  ratio  of  infinitesimals,  is  an  idealized 
concept  that  can  never  be  measured  exactly.  Real  measurements  are 
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Fig,  B.l,  Geometry  of  incident  and  Reflected  Elementary  Beams 
Used  to  Define  the  Bidirectional  Reflectance  Dis¬ 
tribution  Function, 
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always  made  over  some  finite  solid  angle  and  wavelength  interval  and  can 
therefore  yield  only  average  values  Yr  over  those  parameter  intervals. 

The  BRDF  is  basic  in  the  sense  that  all  other  reflectance  or 
scattering  functions  can  be  derived  from  it.  For  example,  Judd  (1967) 
lists  nine  different  kinds  of  reflectance  functions  based  on  the  angular 
extent  of  the  incident  and  reflected  radiation.  All  of  them  can  be 
derived  from  the  BRDF. 

Note  that  the  BRDF  is  a  four- dimensional  quantity  that  can  be 
thought  of  as  an  infinite  family  of  two-dimensional  light  distribution 
functions--one  for  every  possible  angle  at  which  the  incident  beam  can 
strike  the  surface  element.  This  involves  an  overwhelming  quantity  of 
data,  especially  where  high  directional  resolution  is  needed  to 
describe  glints  and  specularities. 


h> 


3 


APPENDIX  C 


I 


-I 


£ 


£ 


I 


COMPUTER  PROGRAM  FOR  INVERSE  SCATTERING  PROBLEM 

The  following  is  a  computer  program  for  calculating  the  effec¬ 
tive  transfer  function,  the  rms  surface  roughness  and  the  surface 
autocovariance  function  from  scattered  light  data.  Once  these  surface 
characteristics  are  known,  the  same  program  can  be  used  to  predict  the 
scattering  properties  at  a  different  wavelength. 
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PROGRAM  SSMA  IN  (INPUT,  OUTPUT,  TAPfc5»  I  NPUT,TAP£M0UTPUT) 

Ini tREw  wonrt 

COMMON  HS(A'-10),  BO(30M),  CSUHNI300),  HUMP),  SCALE,  WAVE,  RMS 
COMMON  oS,  IS,  IA,  ROUTE,  SAMPLE 

DAT*  TYPE  1 ,  TYPES,  TYPE3,  FN0P/4HSIJRF ,  4HAUT0,  4MSCAT ,  4HEN0P/ 

1  OAT*  IFLAU/H/ 

C  *** 

90C*l  FORMAT  (A4,  I  ?  *  A  4  ) 

9990  FORMAT ( I H| , 6 (/ ) , I  0 ( 1 0X , 31 HEWRPR  In  'TYPE'  OF  PROGRAM  DATA,//)) 

C  *  *  »  A********************* *»*»•*«****«*„ ».•»*»*«*. ***,****»*»«»..„»„** 
C  *»*  Rf.  A'l  TYPE  of  I' A T *  tlSFu 

10  RtAi'fS,  90,11 )  SAMPLE,  ROUTE,  TYPE 
IFfTYPE.EO.TvPF  1)  Gil  10  ?U 
IF  tTYPE.E  5.  fYPFS)  GO  TO  3'/* 

IF  f  TYPE  ,EU.  f YP13)  GO  TO  4-1 
IE  fYYPK.tO.EUOP)  GO  TO  09 
WRITE.  (6,949,1) 
r.'J  TO  99 
SO  CAll  SURF 
T F|  *Gt) 

30  Cll  L  AUTO(fFLAG) 

IFI  »Ot/l 
GO  TO  10 
40  f.AU  SCAT 
rO  TO  |0 
99  CitL  EXIT 
Yvl 

SUBROUTINE  SURF 
TlftrpFR  ROUTF 

R{  <L  >,  •'I,  *(30P),  F  ( 4 ) ,  F(5li),  U(i&0) 

COMMON  PS  (SCO),  BOOK!'),  CSUBw  (30K)  ,  H(300),  SCALF,  WAVE,  RMS 

Common  OS,  NTOT,  IA  ,  ROUTE ,  SAMPLE 

EQUIVALENCE  ($,E(1)),  (S1,F(?)),(SS,E(3)),(S3,E(4)) 

DATA  ISTAR/tw*/ 

9'»?0  FORMAT  (A  3,  S,  J*,  13,  lX,FS,l,IX,Et0.9) 

9%’ l  FORMAT (FS . 0 , 1 X , F3 , 1 , 1  A , 13, IF, IS) 

9SM>?  FORMAT  (ShFJ.S) 

99,10  FORMAT ( J  Ml |  ?0{/),  SPY ,  SUBSURFACE  PROFILE  OATA,//, 

I  SOX ,  19‘MITH  WAVELENGTH  OF  ,EiH.4,7H  MICRON,//, 

S  SPY,  SHHOHAPHEO  WITH  A  SCALE  FACTOR  ,  F3, I , //, 

3  SO* ,  1  ‘.-'MS  AMPLE  i"J,,  l  Y ,  A4) 

9901  FflRMATCIH!  ,  1SHSIJRMCF  PROFILE) 

990?  FORMAT!/, Hx,4H«MS»,E9. 3, bH  YmIn»,E9,3,SH  YMAX»,E9,3,//) 

99.*  3  FORfiATtPY,  Ii,  1 « l  A  J  3 

9904  FORMAT  (  /  / , A»X , 9HHIST0GR AM) 

C  **»  HHHHHttlMMttiMHHMHttMUMHtmMtHtOtMOHtllMH 

NTOT*H 

C  A**  OR  IS  SAMPLING  nlSTANCE  ON  PHOTO  IN  MM 
RfcAU (5 * 9(U'0 )  US,  IA,  SCALE,  HAVE 

4  REAL1  (5*9  1C  l )  Ml,  7,  MOfcY,  NOEX 
Z»7/57,S95« 

}>  00  5  1*1,4 

5  Ed)  •  0. 

READ  (5, 900?)  (  A  ( I )  ,  T  *  1  ,  Mt'EX) 

C  ***  PARALLAX  EOUAIIUN 
B»A (1 ) 

C«?»Ml*SINd) 

00  J0  I » I , MOEX 

I 


> 


131 


10  S«54*(I) 

c  ***  DATA  CONVERSION  FOR  SURFACE  PROFILE 
XiMOtK/?. 

V»S/MOEX 
on  aa  i ■ i ,mofx 

S1«>1*(I-X)*(A(I)-Y) 

sa*s?+n-x)*»a. 
aa  S3*.'>!+(A(n-v)'**«» 

M«si/sa 

B«Y-M*X 

do  as  i*i ,moex 

ACI)»A(I)-(M*1*U) 

JOaN'TOT  ♦  I 
as  U(IO)«ACT) 

NTflT  «  Nf'JT  ♦  MOEX 
IF(NOEX,LE.O)  GO  TO  <* 

OS*PS/Ml  a  100(1 

WRJ1E(6,9900)  NAVE,  SCALE#  SAMPLE 
W&m(6,9V0l) 

CALI  PRlNT(i.),OS,NTnT,a) 
c  ***  FIND  TMAX  AND  YMIN 
$a0 

Y1«0(1) 

Ya»c(U 

no  T?  I  *  t , NTOT 
s«s*nn)**a, 

IFCYl.Gf.Od))  GO  TO  3P 

Yi.om 

GO  TO  33 

30  lF(Ya,LE,OU))  GO  TO  35 

vatom 

IS  CONTINUE 

C  «*«  TAH.it  ATE  THF  HISTOGRAM  ANO  CALCULATE  RMS 
X«(Yl-Ya)/N!)EX 
PO  <10  I«t,NOEX 
90  FClJaP,. 

PO  6"  I«l,NrPT 
DO  S«  Kat.NOEX 
L*K-NPFX/a, 

IF(On).GT,L*X)  GO  TO  50 
f  (K)»F  (K)  ♦  ! 

CO  TO  60 
SP  CONTINUE 

tp  continue 

WRITE  (6,99<)<0 

RMS» CS/NTOT) a*, 5 

WRITE  (6,99na)  RMS,  Ya,  Yi 

00  70  I ■ 1 # noex 

NHaf  Cl) 

NRJ  »NH*1 

70  WRI1‘  (6,9903)  NH ,  ( 1  ST  AH  ,<?«*,,  N1*U 
C  .*«*  CALCULATE  THE  SURFACE  AUTOCOVARIANCE 
00  90  I«t,NTOT 
CSUH, ( I ) »Of 
NTT«(NT0T-I)+1 
00  60  K *  1 , NT T 

6P  CSIU’»v.(n«CSUBR(I)  ♦  n(K)*0(K*I-t) 

90  CSUhW(n*CSunw(I)AOS 
RETURN 
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ENO 

SUflROUTINE  AIJTO(IFLAG) 

INTEGER  ROUTE 

Common  HS(3tiO) ,  RQ(3P:»).  CSUR*(300),  H  (300)  ,  SCALE,  WAVE,  RMS 
COMMON  l)S,  IS,  I A  t  ROUTE,  SAMPLE 
DATA  PI/3, U15927/ 

C  *•«  ************** *•**«***•*•••*•*•**•**•*«•*•***«•***•**« ••**••**«•** 

9000  FORMAT(A?,l,l*,EIM,«) 

9001  Fo»MAHM.«,  2(13,  IV),  F2.1,  IX,  Elli.a,  IX,  E9,3) 

90/2  FORMAT  (  8(F9,3,tX)) 

99M0  FORMAT (tMI ,  20(/),  50X,  ?7M$L!FACE  AUTOCOVARIaNCE  OATA,//, 

t  SOX,  19H»ITH  WAVELENGTH  OF  ,E10,a,7H  MICRON,//, 

2  5MX,  ?Hm(.r*phF.I'  wjlil  A  SCALF  F'aCTOR  ,  F3,t,//, 

3  sox,  10HSAHPIF  NO.,  IX,  A«) 

9901  F0RM*T(1H1  ,  IRHSCATTf.OlNG  FUNCTION) 

99H2  FORMAT  (//,?9M  IMF.  PRS  SIRFACF  ROUGHNESS  IS,  F9.3,  7H  MICRON) 

9903  FOR’IAT(jHi  ,  <i2H(SVSTph)  FOUR ti-R-RESSFL  TRANSFORM  FUNCTION) 

990H  FORMAT ( t H 1 ,  3 1 HSUW FACE  AUTOCOVARIANCE  FUNCTION) 

990S  FORMATUHl  ,  17HTRANSFER  FUNCTION) 

C  ***  **«**.. •***.•«.•*•*.*.•.**«.•.«**••.*..*.*****•*.«*. ************** 
IF  (ROUTE  O)  GO  TO  <1 
RtAO(S,900O)  SCALE,  WAVE 
WR I TF (b, 9900)  nave,  SCALE,  SAMPLE 
8  IFUrtAG.eO.l)  GO  TO  5 

RE AO (S, 9001 )  OS,  IS,  IA,  SCALE,  WAVF,  RMS 
RE  AT  (5, 9.1 1? )  (rSUBW(t),  1*1, IS) 

WRIU(*»,99'V)  WAVE,  SCALE,  SAMPLE 
S  »RI 1 E (h ,  99041 ) 

WRITE  (b,99l*2)  RMS 

CALL  PRIN? (CSUHW,03» IS, 2) 

ns*OS/wAVE 

C  ••*  CALCULATE  THE  TRANSFER  FUNCTION 
00  10  1*1, IS 

10  H(T)«tXP(  (9 *P  l /WAVE )**2.*tCSUiml)-RMS*»2.)> 

WRI TC (b , 9905) 

CALI  PRTNT(h,DS,IS,1) 

C  *«*  CALCULATE  f,  FROM  RMS  SURFACE  ROUGHNESS 
B«1  ,  •  1  ,  /  E  X  P  (  (R<*S*R»PI  / WAV  (•')**?,) 

C  *«*  F I  NO  THE  (SYSTtM)  FOUR ItR-St SGEL  TRANSFORM  FUNCTION 
00  20  1*1, IS 

20  Bil(l)*H(I)*(R«i,) 

WR I  TC (b , 9903) 

CALL  PRINT (HU,0S,IS,1) 
c  **•  OETEPHTNF.  THE  OUTPUT  JNTFPVAL  OA 
DA*, 02 

C  ***  FINO  THE  SCATTERING  FUNCTION 

CALL  HANKf.L(R'),OS,  IS,B3,0A,IA) 

WRI1F(6, 99-51) 

CAll  PRINT  (MS,  t'A  ,  I  A,  3) 

c  ***  F I  NO  THE  (SYS  f F  M)  FOURIEP-FESS' L  TRANSFORM  FUNCTION 
XM*IS-1 
00  3P  T •  1 , 1 S 
X»I-1 

30  B0(I)*H0(I)*EXP(«PT*(X/XM)**2,) 

WR  fl F (b,9903) 

CAU  PRIHT(tlO,OS,IS,l) 

C  ***  FIND  THE  SCATTERING  FUNCTION 

CALI  HANKEL(MU,0S,I5,B3,t)A,  I  A) 

WRITE (b,990l ) 

CALL  PRINT (US , OA , IA,3) 


return 

eno 

SUBROUTINE  SCAT 
INTEGER  ROUTE 

COMMON  B$(3UO),  BQ(3&0),  C»UBw(30B),  M(3BtO,  SCALE,  WAVE,  RMS 
COMMON  ITS,  IS,  IA,  ROUTE,  sample 
DATA  Pl/3, 1415927/ 

C  aaa  aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa aaaaaaa aaaaa aaaaaaaaaaaaaaaaaaaa 

90UI  FUMMIF},?,  IX,  2(13,  IX),  F?,l,  IX,  E1P.4) 

9PU?  F()RM4T(  P(F9.3,1X)) 

9903  FORMAT (1M1  ,  3K(/),  5PX,  24HSCAT  TER  INC  FUNCTION  DATA,//, 

1  SOX,  19MWITM  WAVELENGTH  OF  ,K1P.4,7H  MICRON,//, 

2  SOX,  2SHGKAPHF.I'  WITH  A  SCALE  FACTOR  ,  F3, I , //, 

3  bOX ,  1CHSAMPLE  NO,,  IX,  A4) 

9901  FORMAT ( 1  H  J , 19HSCATTFRIM,  FUNCTION) 

9903  FORMAT(//,?«h  ThE  Rms  SURFACE  pnijGHNFSS  IS,  F.9,3,  7N  MICRON) 

9903  FORMAT ( t H | ,  9?M (SYSTF  M)  FOUR 1FR-HF  SSEL  TRANSFORM  FUNCTION) 

99F4  FORMAT  ( 1 M 1 ,  3 1 HSURF ACE  AU T Otnv AK 1 AMCE  FUNCTION) 

99Fb  FORMAT (IMl,  PshSYSTEM  TRANSFER  FUNCTION) 

RE  AO  ( «j ,  90v!  1 )  OA,  IA,  IS,  scale,  WAVE 
READ  (5,  4t'?>?)  (MS  (  I  )  ,  1*1,  I  A  ) 

WRITF  (S,990t:)  RAVE,  SCALE,  SAMPLE 

WRITE (F, 99011 

CALI  PRINT (MS,OA, IA,3) 

f.  aaa  ORTf. Rm I NE  Thf  OUTPUT  INTERVAL  OS  (ASSUMING  AR  IS  NON-INCREASING) 
CALL  INTFR(wS,OA, IA,h$,lS,0) 

IF (ROUTE, EO.O)  0S«.2/AAVF 
OSiORaSCALE 

c  aaa  FINO  THE  (SYSTEM)  FOIIRIER-OESSEL  TRANSFORM  FUNCTION 
CALL  HAN*EL(PS,I)A,IA,B0,0S,I3) 
wRITF (*>,990.3) 

CALL  PRlNT(HO,l)S,IS,t)- 
C  AAA  C  At  CliL  ATE  THr  TRANSFER  FUNCTION 
BtROf  1) 

00  20  I • t  ,  IS 
20  M(I)«(!,-B)  ♦  H Q ( I ) 

WRITF (A, 9905) 

CALL  PRINT (H, OS’,  IS,  1 ) 

C  AAA  CALCUIATF  SJGmA  SNUAHEO  ANO  RMS 

SIGS0*(wAVF/(4,aFI))aa2,aAL0G(1,/(1 ,-B)) 

RMS  «  SIGSQaa  , b 

C  aaa  CALCULATE  THF  SURFACE  AIJ TOCOVaRI ANCE  FUNCTION 
00  30  T»t,!S 

3/  CSUbt  (J)i  (wAVf /(«,aPI))**2»*AL0G(H(I))  ♦  SIGSO 
L'S  •  DSaWAVE 
WRITF (A, 9904) 

WRITF (4,9902)  RMS 

CALL  PRINT (CSU»w,OS, IS, 2) 

IFCRPUTE.60,1)  CALL  AUTOU) 

RETURN 

EnO 

SUBROUT  INF  INTt  R(X,CX,IX,OZ,a,IFLAG) 

REAL  X(3P,*.),  075(5) 

DATA  EPS/ ,01 /,  07$/, 1 »  ,2,  ,?5,  ,5,  »,/ 

I M  A  X  *  ? 

VMAX«X(1) 

00  b  I  *2 , t  X 

TF(TMAX,GF.X(l))  go  to  s 
VMAX1XU) 


X 

> 


IMAX*1 
5  CONTINUE 
l)Z»YPAX*EPS 

If  (If LAG.t-Q.O  on  TO  ?5 
DO  10  I*TMAX,IX 

f  IKxm.lT.OZ)  GO  TO  20 

tit  CONTINUE 

?.1  OZ*l,/(I-n/nx/I2«3, 

iK!a*(,,tcj,a)  oz»t,/(i-!)/ox/i2/3, 

25  B«  AI.OGlKOZ) 

K»  I N T  ( o T  *•  I 

IKP.LT.P)  K*INT(B) 

TES1 *  10.**  CB-K) 

|  on  3?  1*1,5 

IKTFST.Lt.OZSU))  GO  TO 
30  CUNT/NIIE 
<tn  02*  nzsci)*i0.**K 

RETURN 

ENO 

SUBROUTINE  PRInT(h,0X,  IH,IU) 

REAL  P(3VM),  PTttOl),  UN  ITS  (3) 

$  OA 1  A  UNI  fb/ftHI.  AMHliAt  MU*  IDRON,hH  BETA  / 

c  *»«  *•***«*****»•**«»*»»•*»»*. *..*..«».**,*.«****'***.**» 

9001  F0Rl'ATf//,2M,«M».l'  ,B(12X,EH,2)  /  «X  ,  2HX  (  ,  A6,  1  H ) ,  3X ,  9HAMPL  I TUOE . 

C  UX»lMA#5(il5n»*«“>)»5  »•-•«»*)) 

99.12  FORMA  1  («X,E9,3,3X,F9.3.«X,  Id)  AO 

C  ***  ***«*»*. A********************.*********************** 

no  in  1*1,101 
to  PT ( I } *  1 H 

y  X*f! 

C*l  L  1NTFR(W,0X*  I W ,  DAMP ,  |  ?i(l ,  1  ) 

WRITE  (f.,99('l )  (OAMPAl,  I*?P, 100,20),  UNtTS(IU) 

00  Jei  l*t,lw 
PT(1).1H( 

PT  (  1511  )  *1  111 

K  »A0 A ( w ( I ) /DAMP)  ♦  t,S 

v  IKK. GT. ini)  K*121 

l  PT(K).tH» 

IKWU).LT.(I)  PT  (K )  » 1  HO 
WKJU (6,9902)  X,*(I),PT 
PT (K) • *  H 

x*x*ox 

30  CONTINUE 
RETURN 

v  End 

*  5UBP0UTIN1  l'ANktKW,OR,lW,V,r>kMO,IV) 

REAL  Wfiji.o,  VO  no) 

OAT*  C/E.?8318'»«/ 

00  an  l«i, IV 
A«0, 

RHO* ( 1 • 1 ) *ORHO 
00  IB  I»*2,  lw 

v  R«0R*(IR-2)A0«/a. 

*  in  A»A*R»('«,(JRMKIR-l))*MJl>(C«RHQ*R) 

20  VU)*A»n«»C/2, 

Rf  TURN 
END 

FUNCTION  SJB(X) 

IF(X,LT,0,)STOP 
BJO*  1  ,0 
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IF(X.EO.{!,)»ETURN 
CONST*?. /X 
M  ■  t  tJ  ♦  1 ,3»X 
pMt«i,«F-ae 
fm«b. 

AUF«H. 

HI »M-1 
M2»M«2 

oo  ao  K*t,H?,a 

BMK«  »Fnl  *CONST-FH 

FM«FMl 

FMJ*6MK 

8mK*(m1-)<)*FM1*C0NST-FH 
FMif  Ml 
FM1 «HMK 

?0  »LF«*LF»0HK 
AUF«?,0»*UF«flMK 

bjo«»mk/*uf 

RETURN 

EN0 


APPENDIX  D 


COMPUTER  PROGRAM  FOR  PROCESSING  SURFACE 
PROFILE  DATA  FROM  ELECTRON-MICROGRAPH  STEREOCOMPARATOR 

The  following  is  a  computer  program  for  processing  surfoce 
profile  data  from  the  electron-micrograph  ste/eocomparator.  It 
provides  the  surface  height  distribution,  the  rms  surface  roughness. 


and  the  surface  autocovariance  function. 


STFR3 


13:06PDT 
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10  DIM  A(200).  F<200) 

20  RFAO  M#  Z.  N>  HI 
30  Z=Z/S7.29S8 
AO  FOR  I  =  1  TO  N 
SO  RFAO  A  (  I  ) 

60  NFXT  1 
70  R  =  AM) 

80  0  =  2*M*SIN(Z) 

90  FOR  1  =  1  TO  N 

ioo  Am  =  c a ( i )  -  p)/r 

1 10  NFXT  I 

120  S  *  SI  =  S2  =  S3  =  0 
130  FOR  I  =  I  TO  N 

1  AO  S  =  S  +  A  (  I  ) 

ISO  NFXT  I 

160  Y  »  S/N 

170  X  *  N/? 

180  FOR  I  =  1  TO  N 

190  SI  s  SI  ♦  (I  -  X)*(A(I)  -Y) 

200  SP  =  SP  ♦  (I  -X) iP 

2  10  S3  =  S 7  ♦  <A(  I  )  -  Y) »P 
2P0  NFXT  I 

230  M  *  SI  /S2 

2A0  P  *  Y  -  H*X 

250  C  =  S1/S0R<S?*S3> 

260  S  «  Y!  =  YP  =  0 

270  FOR  I  =  1  TO  N 

280  A(I)  =  A ( I )  -  (N* I  +  B > 

290  S  *  S  ♦  A(I) *2 

300  IF  ri  >  A(I5  THFN  330 

310  Y1  =  A  {  I  ) 

320  GO  TO  3S0 

330  IF  Y2  <  A  f  I  )  THFN  3S0 

340  YP  =  A ( l ) 

3S0  NFXT  I 

360  PRINT  "RPS" #  "HFRIT  FUNCTION".  "YNAX". 
370  PRINT 

380  PRINT  SORf S/N) .  C.  Y1.  Y2 
390  0  =  (Y1  -  Y2)/K1 
AOO  FOR  I  =  1  TO  K1 
A10  F(I)  =  0 
A20  NFXT  I 


"YMIN" 
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430  FOR  I  =  I  T0  N 

440  FOR  K  =  1  TO  Kt 

450  J  -  K  -K»/2 

460  IF  A(I)  >  J*D  THEN  490 

470  FCK)  =  F  OO  ♦  1 

480  GO  TO  500 

490  NEXT  K 

500  NFXT  I 

510  PRINT 

520  PRINT 

530  PRINT  "  HF1GHT  DISTRIBUTION" 

540  PRINT 

550  FOR  I  =  1  TO  K1 
560  PRINT  F( I ) 

570  F(I)  =  0 
580  NFXT  I 
590  PRINT 
600  PRINT 

610  PRINT  "N" «  "Y",  “AIJTOOOVAR 1 ANCF" 
620  PRINT 

630  FOR  I  =  I  TO  N 
640  FOR  K  =  1  TO  CN-1 > 

650  F < I )  =  F  <  I  )  ♦  A  ( K  )  *  A  <  K  ♦  1  -  I) 
660  NFXT  K 

670  PRINT  I.  4(1).  F( I ) 

680  NEXT  I 


690 

DATA 

toooo 

.  4.  40,  11 

♦ 

700 

DATA 

2.^0. 

2.87. 

2.93. 

3.05. 

3.05. 

3.12. 

3.12. 

3.13, 

3.25 

710 

DATA 

2.90. 

2.67. 

2.80. 

2.83.3 

.00,  2 

.90.  2 

•  23.  9 

•  87,  2 

.63 

720 

DATA 

2.86. 

2.65. 

2-76. 

2  -72. 

2.78. 

2.78. 

2.96, 

2.82. 

2.93 

730 

DATA 

2.81  . 

3.00. 

2*90. 

2.83. 

2.89, 

2.90. 

2.98, 

2.86, 

3.05 

7  40 

DATA 

2.88. 

3  .03. 

2.88.3 

.29 

750 

DATA 

2.47, 

2.33. 

2.34, 

2  .38. 

2.28. 

2.33. 

9.43. 

2.48. 

9.61 

760 

DATA 

2.58. 

2.53. 

2.35. 

2.36. 

2.30. 

2.20. 

2.15, 

2.22. 

9.16 

770 

DATA 

1.84, 

1  .92. 

2.09. 

2  .93. 

2*27, 

2.31, 

9.19. 

1  .95. 

2.05 

780 

DATA 

2.00. 

1  .94, 

2*10. 

2.23. 

9.21, 

2.00. 

1  .89. 

9.12. 

2.90 

790 

DATA 

9  .03. 

2.00. 

2.03. 

1  .96 

800 

DATA 

2.87. 

2.84, 

2.85. 

2.71  . 

9.94, 

2.87, 

3.08. 

3.10. 

3.95 

810 

DATA 

3  .00. 

2.97. 

3.14. 

2.98. 

3.00. 

2.98. 

3  .06, 

2.70. 

9.87 

820 

DATA 

2.80. 

2.74, 

2.79. 

3.00. 

3.20. 

2.97, 

3.11, 

9.76, 

9.86 

830 

DATA 

2.74, 

9.83. 

3.01  . 

3.16. 

3.06. 

2.79, 

9.73, 

9.85. 

2.90 

840 

DATA 

2.73. 

2  .93, 

2.90. 

3  .00 

850 

DATA 

2.00. 

2.08. 

2.93. 

2  .02. 

2.22. 

2.20. 

2.14, 

2.24, 

2.30 

860 

DATA 

2.30. 

2.39. 

2.31  . 

2.21 . 

2.42, 

2.38. 

2.38. 

2.37. 

9.33 

870 

DATA 

9,29. 

2.27. 

2.01  . 

2.19.9 

.28.  2 

.27,  2 

.21.  2 

.34,  2 

.08. 

880 

DATA 

2.25. 

2.27, 

2.15. 

2.34, 

2.40, 

2.22. 

2.15 

.  2.16 

,  2.19 

890 

DATA 

2.15. 

2.31  . 

2  •  44  * 

2.52 

999  FNO 
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DESCRIPTION  OF  SCANNING  FECO  INTERFEROMETER 
USED  FOR  DETERMINING  STATISTICAL 
PROPERTIES  OF  OPTICAL  SURFACES 


The  following  is  the  manuscript  of  a  paper  presented  at  the 
1974  Annual  Meeting  of  the  Optical  Society  of  America  in  Houston, 
Texas  (J.  Bennett,  1974). 
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Use  of  Interferometry  for  Determining  the  rms  Roughness,  Autocovariance 
Function,  and  Other  Statistical  Properties  of  Optical  Surfaces* 

Jean  M.  Bennett 

Michelson  Laboratory,  Naval  Weapons  Center 
China  Lake,  California  93555 

ABSTRACT 

A  FECO  Scanning  Interferometer  will  be  described  that  can 
measure  very  small  height  differences  with  a  lateral  resolution 
of  2  microns  to  yield  statistics  for  optical  surfaces. 


The  other  papers  in  this  Symposium  on  Techniques  in  Surface 
Interferometry  have  mainly  considered  the  use  of  interferometers 
to  study  the  contours  on  optical  surfaces,  or  more  specifically 
the  deviations  from  the  desired  surface  contours.  These  devia¬ 
tions  are  sometimes  called  the  figure  of  the  surface,  and  opti¬ 
cians  frequently  talk  of  half-wave  or  quarter-wave  optical 
surfaces  when  they  mean  that  the  deviations  from  a  perfect  plane 
or  curved  surface  are  one-half  or  one-quarter  of  the  wavelength 
of  the  light  used  to  test  the  surface  (traditionally  the  mercury 
green  line  at  5461  X) .  The  figure  of  an  optical  surface  is 
important  in  determining  the  resolving  power,  focusing  proper¬ 
ties,  and  aberrations  in  an  optical  system.  For  this  reason 
much  work  has  been  devoted  to  interfacing  the  interferometer, 
which  can  sense  figure  errors,  to  the  polishing  machine  which 
can  eliminate  them. 

Optical  technology  has  now  progressed  to  the  point  where 
another  parameter,  the  microroughness,  also  becomes  important. 
Microroughness  on  the  surface  scatters  some  of  the  light  into 
unwanted  directions.  Scattering  cannot  only  reduce  the  contrast 
in  optical  images  by  removing  light  from  the  bright  areas  and 
filling  in  the  dark  ones,  but  it  can  also  drastically  reduce 
the  optical  throughput  of  a  system  long  before  the  resolving 
power  is  affected.  Scattering  is  also  a  serious  problem  when 
one  is  trying  to  observe  a  weak  object  that  is  very  close  to 
a  bright  object.  This  situation  occurs  frequently  in  astronomy 
when,  for  example,  one  is  trying  to  observe  details  in  the  solar 
corona  near  an  occulting  disk,  or  when  looking  at  a  faint  star 
located  close  to  the  moon  or  a  bright  star. 

The  problem  of  scattering  from  optical  surfaces  has  been 
around  for  a  long  time.  Even  before  1900  Albert  A.  Michelson, 
the  man  best  remembered  for  his  measurements  of  the  velocity  of 


^Reproduced  here  with  permission  of  the  author. 
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light  and  the  Michelson-Morley  experiment,  thought  about  light 
scattering.  At  the  Michelson  Museum  in  China  Lake  we  found 
an  entry  scribbled  in  cne  of  Michelson's  pocket  notebooks 
along  with  shopping  lists,  a  prescription,  notes  for  a  coming 
lecture,  possible  causes  of  error  in  an  experiment,  and  other 
diverse  items.  The  entry  is  shown  in  Slide  1  and  reads,  "Find 
relation  between  roughness  of  surface  and  angle  of  scattering. 
Much  has  happened  in  the  80  or  so  intervening  years,  but  we 
still  have  not  satisfactorily  solved  the  problem  Michelson 
hastily  noted  down  around  1890. 

What  we  have  learned  about  the  relation  between  surface 
roughness  and  scattered  light  is  that  the  total  hemispherical 
scatter  from  a  surface  (i.e.,  all  the  light  scattered  into  a 
hemisphere)  is  related  primarily  to  the  heights  of  the  surface 
irregularities  when  these  heights  are  small  compared  to  the 
wavelength  of  light.  By  measuring  the  total  hemispherical 
scatter  (frequently  called  TIS  for  short)  and  assuming  a 
Gaussian  distribution  for  the  heights  of  the  irregularities, 
we  can  obtain  a  value  for  the  rms  roughness  of  the  surface. 

For  many  types  of  surfaces,  particularly  polished  glass,  fused 
quartz,  calcium  and  magnesium  fluorides,  etc.,  this  rms  rough¬ 
ness  value  is  in  excellent  agreement  with  the  roughness  value 
obtained  from  interferometric  measurements,  a  technique  I  will 
describe  in  detail  in  a  few  minutes.  However,  where  the  scat¬ 
tering  theory  falls  short  is  in  predicting  the  effects  of 
scattered  light  from  unusual  optical  surfaces  such  as  polished 
alkali  halides,  polished  metals,  electropolished  metals,  and 
micromachined  metals.  In  these  cases  TIS  measurements  can 
yield  an  effective  rms  roughness  value  that  is  either  consid¬ 
erably  smaller  or  considerably  larger  than  the  interferometri- 
cally  measured  value.  More  troubling  even  than  this  is  the 
observation  that  we  are  not  able  to  correctly  predict  the 
angular  dependence  of  scattered  light  about  the  specular  dir¬ 
ection  even  for  the  smoothest  polished  glass  surfaces.  This 
situation  arises  because  all  the  scattering  theories  assume 
a  Gaussian  autocovariance  or  autocorrelation  function  for  the 
surface  and  none  of  the  real  surfaces  we  have  studied  have 
Gaussian  autocovariance  functions.  For  this  reason,  in  order 
to  have  a  theory  which  correctly  predicts  the  effects  of  scat¬ 
tering  from  a  surface,  we  need  to  know  the  autocovariance  func¬ 
tion  and  other  statistics  of  the  actual  surface.  How  to 
measure  these  statistical  properties  is  the  subject  I  am  going 
to  discuss  for  the  remainder  of  my  talk. 

Interferometry  has  been  shown  to  be  an  excellent  method 
for  looking  at  very  small  height  differences  on  surfaces,  height 
differences  of  a  few  angstroms,  i.e.,  a  few  thousandths  of  the 
wavelength  of  light.  Tolansky  pioneered  this  type  of  inter¬ 
ferometric  technique  and  gave  the  interference  fringes  the 


enigmatic  name  of  Fringes  of  Equal  Chromatic  Order,  or  FECO  for 
short.  The  two  main  things  that  are  important  about  this  type 
of  intereference  fringes  are  (1)  that  they  contour  height  varia¬ 
tions  on  surfaces,  and  (2)  they  occur  as  wiggly  black  lines  in 
a  continuous  spectrum  of  light  reflected  from  the  interferometer. 

1  am  making  the  second  point  to  distinguish  FECO  fringes  from 
the  more  common  Fizeau  fringes  which  also  contour  irregularities 
on  optical  surfaces,  but  which  are  formed  in  monochromatic  light. 
FECO  fringes  have  advantages  over  Fizeau  fringes  in  that  the 
order  of  interference  of  the  fringe  is  always  known,  and  small 
areas  of  the  surfaces  can  be  studied  at  will  without  having  to 
readjust  the  tilt  of  the  interferometer  plates. 

I  am  now  going  to  describe  a  FECO  system  we  have  built  to 
measure  the  statistical  properties  of  various  types  of  optical 
surfaces  and  I  will  show  you  samples  of  the  data  we  have  obtained 
with  this  system.  I  will  also  mention  some  types  of  experiments 
we  are  planning  for  the  future. 

A  photograph  of  the  FECO  Scanning  Interferometer  is  shown 
in  Slide  2.  Most  of  the  instrument  consists  of  a  signal  averager, 
minicomputer  and  teletype,  and  the  interferometer,  the  heart  of 
the  experiment,  is  the  smallest  part.  The  optical  arrangement 
is  shown  in  Slide  3,  The  interferometer  I  consists  of  the  sample 
to  be  studied,  coated  with  an  opaque  layer  of  silver  (upper 
plate)  and  a  super  smooth  surface  of  polished  fused  quartz 
coated  with  a  semi-transparent  film  of  silver  of  approximately 
95%  reflectance.  The  two  optical  surfaces  are  very  close 
together,  being  separated  by  only  a  few  half  wavelengths  of 
light.  The  actual  spacer  consists  of  the  dust  particles  on 
the  two  surfaces.  The  interferometer  is  illuminated  in  reflec¬ 
tion  by  a  collimated  beam  of  white  light  from  a  xenon  arc.  The 
important  feature  of  the  FECO  system  is  lens  L2  which  focuses 
an  image  of  the  interferometer  surfaces  on  the  slit  S  of  a  con¬ 
stant  deviation  spectrograph.  Thus,  the  interference  fringes 
which  are  viewed  in  the  focal  plane  of  the  spectrograph  contour 
the  irregularities  on  the  pair  of  optical  surfaces.  A  picture 
of  what  might  be  observed  is  shown  in  the  circular  inset. 

Three  mercury  lines  are  included  for  wavelength  calibration 
purposes,  but  the  information  about  the  surface  topography  is 
given  in  the  wiggly  interference  fringe.  There  is  a  one-to- 
one  correspondence  between  the  wiggles  on  the  interference 
fringe  and  height  variations  on  the  pair  of  optical  surfaces. 

To  get  an  idea  of  the  magnitudes  of  the  quantities  we  are  deal¬ 
ing  with,  the  wavelength  variations  can  yield  information  about 
variations  of  the  heights  of  irregularities  of  the  order  of  a 
few  angstroms.  The  lateral  resolution  is  much  smaller,  so  that 
the  length  of  the  interference  fringe  corresponds  to  a  distance 
of  one  mm  on  the  interferometer  surface.  The  width  of  the 
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spectrograph  slit  determines  the  other  dimension,  so  we  are 
actually  looking  at  an  area  one  mm  long  by  about  2  microns 
wide. 


To  obtain  the  statistics  for  the  section  of  the  inter¬ 
ferometer  surface  we  are  observing,  we  scan  the  spectrum  line 
by  line  using  a  special  slow  scan  TV  camera.  After  suitable 
ageraging,  the  information  on  one  scan  line  appears  as  shown 
in  the  upper  right  hand  part  of  the  slide.  The  desired  infor¬ 
mation  is  the  wavelength  of  the  portion  of  the  interference 
fringe  included  in  the  scan  line  and  this  is  obtained  from 
the  computer  analysis.  The  information  from  the  entire  frame 
consists  of  wavelengths  of  512  equally  spaced  points  on  the 
interference  fringe,  so  we  have  data  for  512  different  areas 
on  the  surface,  each  one  of  which  is  a  little  square  2  microns 
on  a  side.  In  the  statistical  analysis  a  least  squares 
quadratic  curve  is  calculated  from  all  512  wavelengths  and 
defines  the  mean  surface  level.  Then  wavelength  differences 
from  this  curve  are  converted  into  height  differences  above 
and  below  the  mean  surface  level.  Using  the  height  differences 
we  can  determine  the  rms  roughness,  height  distribution  func¬ 
tion,  slope  distribution  function,  rms  slope,  autocovariance 
function,  and  other  statistical  parameters  for  the  surface. 

I  am  now  going  to  show  some  data  that  are  typical  of  what 
we  have  obtained  for  various  types  of  surfaces.  We  have 
studied  very  smooth  glass-type  surfaces  such  as  fused  quartz, 
Cervit,  calcium  fluoride  and  magnesium  fluoride,  polished 
alkali  halides  (potassium  chloride  and  sodium  chloride), 
polished  metals  (copper,  beryllium  copper,  titanium  and  titan¬ 
ium  alloys,  molybdenum,  and  stainless  steei)>  electropolished 
nickel,  machined  copper,  and  holographic  gratings.  All  of  the 
smoothest  polished  glass-type  materials  have  similar  statistics, 
and  Slide  4  shows  results  for  one  of  these.  This  is  an 
extremely  smooth  calcium  fluoride  surface  polished  by  Abe 
Klugman  of  the  Northrop  Corporation  and  had  a  visually  measured 
roughness  of  9.4  X  rms.  At  the  top  of  the  slide  is  a  Polaroid 
photograph  of  the  interference  fringe  and  directly  below  it  is 
a  TV  scan  of  the  center  line  of  the  fringe.  Note  that  the 
wiggles  on  the  fringe  represent  height  differences  of  consid¬ 
erably  less  than  10  X.  The  autocovariance  function  shown  below 
can  be  roughly  considered  as  the  correlation  between  points 
on  the  surface  separated  by  the  amount  shown  on  the  x  axis. 

For  the  very  smooth  surfaces  there  is  positive  correlation 
between  closely  spaced  points  but  those  farther  away  are 
random. 

The  height  and  slope  distribution  functions  for  the  same 
surface  are  shown  on  Slide  5.  Note  that  both  measured  distri¬ 
bution  functions  (the  histograms)  are  very  close  to  Gaussian, 
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and  actually  the  smooth  curves  are  Gaussians  having  the  same 
areas  under  the  curves  as  do  the  measured  ones.  The  slope 
distribution  function  is  only  half  a  Gaussian  curve  because 
we  do  not  distinguish  between  positive  and  negative  slopes. 

Polished  alkali  halide  surfaces  and  some  polished  metal 
surfaces  are  similar  in  that  both  are  composed  of  macro  and 
microscratches  with  no  smooth  areas  in  between.  In  Slide  6 
we  see  the  scanning  camera  trace  for  a  polished  KC1  surface 
and  the  autocovariance  function.  Note  the  oscillations  in 
the  autocovariance  function  which  indicate  longer  range  cor¬ 
relations  than  those  observed  for  the  smoothest  surfaces. 

Note  also  that  neither  this  autocovariance  function  nor  the 
preceding  one  were  Gaussian  in  shape.  In  fact,  we  have  never 
observed  an  autocovariance  function  that  did  have  a  Gaussian 
shape.  In  Slide  7  we  see  that  there  is  slight  asymmetry  in 
the  height  distribution  function  although  the  slope  distribu¬ 
tion  function  seems  to  be  a  very  good  Gaussian.  Incidentally, 

I  should  point  out  that  the  rms  roughness  values  shown  on  the 
height  distribution  function  histograms  are  about  a  factor  of 
two  smaller  than  the  visually  measured  values.  This  is  because 
the  scanning  camera  does  not  take  into  account  the  width  of 
the  fringe,  only  its  center  line.  Visually  we  measure  the 
extreme  width  of  the  fringe  and  convert  the  peak- to-val ley 
roughness  to  an  rms  value  by  dividing  by  2 f2.  There  is  a  very 
good  linear  relationship  between  visually  measured  roughnesses 
and  the  values  obtained  from  the  scanning  camera. 

In  Slide  8  we  see  data  for  molybdenum,  a  typical  polished 
metal  surface.  The  autocovariance  function  with  the  oscilla¬ 
tions  is  similar  to  that  for  KC1.  The  height  and  slope  distri¬ 
bution  functions  shown  in  Slide  9  are  reasonably  good  Gaussians 
although  there  is  some  raggedness  on  the  height  distribution 
function. 

Recently  very  low  scatter  electroless  nickel  mirrors  have 
become  available.  In  Slide  10  is  shown  the  autocovariance 
function  for  one  of  these,  which  had  a  visually  measured  rough¬ 
ness  of  23.6  X  rms.  This  surface  is  a  gradually  undulating 
one  with  almost  no  obvious  scratches.  The  autocovariance  func¬ 
tion  also  has  a  lower  frequency  oscillation  than  was  observed 
for  the  KC1  and  molybdenum  surfaces.  In  .Slide  11  we  see  that 
there  is  a  definite  asymmetry  in  the  height  distribution  func¬ 
tion.  There  are  proportionately  more  small  bumps  on  the  sur¬ 
face  than  there  are  small  holes.  However  the  slope  distribu¬ 
tion  function  is  an  extremely  good  Gaussian. 

One  of  the  most  interesting  surfaces  we  have  encountered 
are  those  made  on  a  special  type  of  lathe  by  a  single  point 
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diamond  tool.  Recently  we  determined  the  statistics  on  one 
which  had  a  visual  roughness  value  of  only  22.5  X  rms. 

Slide  12  shows  that  the  autocovariance  function  for  this  sur¬ 
face  has  only  a  long  period  oscillation.  We  probably  are  not 
resolving  the  individual  grooves  made  by  the  diamond,  but  only 
multiples  thereof.  The  height  and  slope  distribution  func¬ 
tions  on  Slide  13  are  reasonably  good  Gaussians  even  though 
there  must  be  obvious  periodicity  in  the  surface. 

As  a  final  set  of  statistics,  we  see  in  Slide  14  a 
tracing  of  a  holographic  grating  with  a  nominally  sinusoidal 
groove  shape  formed  in  a  photographic  emulsion.  Tnis  sample 
was  kindly  furnished  to  us  by  John  Stover  of  Dow  Chemical 
Company,  Rocky  Flats,  Colorado.  The  autocovariance  function 
is  for  only  one  scan  and  clearly  shows  the  periodicity  of  the 
surface.  The  height  distribution  function  on  Slide  15  is 
clearly  not  Gaussian,  and  I  did  not  even  attempt  to  put  a 
Gaussian  curve  through  the  data. 

Our  plans  for  the  future  include  increasing  the  sensi¬ 
tivity  of  the  scanning  camera  so  we  can  obtain  scans  from  many 
more  portions  of  the  surface.  Most  of  the  data  shown  here  are 
averages  of  8-10  separate  scans.  We  also  hope  to  automate 
the  scanning  interferometer  so  that  we  can  make  equally  spaced 
scans  adjacent  to  each  other  to  obtain  statistical  data  on 
a  square  one  mm  on  a  side. 

In  conclusion,  I  have  described  a  type  of  interferometer 
which  can  be  used  to  obtain  statistical  information  about  the 
topography  of  optical  surfaces.  This  instrument  can  distin¬ 
guish  height  differences  of  only  a  few  angstroms  and  has  a 
lateral  resolution  of  about  2  microns.  Using  the  statistics 
of  the  surface,  we  hope  to  be  able  to  predict  the  scattering 
properties  of  the  surface. 


fcrproductlo©  of  oo  original  itmu 
Id  tba  collection  of 
DOROTHY  MICKKL3GK  LIYDW6T0* 


Tht  firat  of  two  anall 
pocket  notebook*  of 
Klchelaoo  Memoranda 


SLOPE 

DISTRIBUTION  FUNCTION 


ELECTROLEU  NICKEL 


SLOPE  (X /micro n) 


SLIDE  11 


MACHINED  COPPER 


MACHINED  COPPER 


HEIGHT 

DISTRIBUTION  FUNCTION 


SLOPE 

DISTRIBUTION  FUNCTION 


RMS 

ROUGHNESS 


DISTANCE  FROM  MEAN  SURFACE  LEVEL  (A) 

SLOPE  (A/fmeron) 

SLIDE  13 

holographic  grating 

l\v  SPACING 

autocovariance 

function 


0.2 

LENGTH  p(mm) 


SLIDE  14 


HOLOGRAPHIC  GRATING 
21*/  SPACING 


HEIGHT 

DISTRIBUTION  FUNCTION 


DISTANCE  FROM  MEAN  SURFACE  LEVEL  (I) 


SLOPE  (X/mleron) 


SLIDE  15 


LIST  OF  REFERENCES 


Ament,  W.  S.,  "Toward  a  Theory  of  Reflection  by  a  Rough  Surface," 

Proc.  IRE  41:142  (1953) . 

Ament,  W.  S.,  "Reciprocity  and  Scattering  by  Certain  Rough  Surfaces," 
IRE  Trans.  AP-8: 167-174  (1960). 

Barrick,  D.  E.,  "Two  Experiments  Yielding  Lunar  Surface  Information 
Employing  Polarized  Radar  Waves,"  Ohio  State  University  Research 
Foundation,  Columbus,  Ohio,  Available  from  the  National  Technical 
Information  Service,  No.  N66-14962  (1965). 

Barrick,  D.  E.,  "Rough  Surfaces,"  Radar  Cross  Section  Handbook, 

Vol.  2,  G.  T.  Ruck,  et  al.,  eds..  Plenum  Press,  N.  Y.,  pp.  671-772 
(1970). 

Bass,  F.  G.,  "Very  High  Frequency  Radiowave  Scattering  by  a  Disturbed 
Sea  Surface,"  IEEE  Trans.  AP- 16: 554-568  (1968). 

Beaglehole,  D.,  "Optical  Excitation  of  Surface  Plasmons,"  IEEE  Trans. 
ED- 17: 240-244  (1970). 

Beaglehole,  D.,  and  0.  Hunderi,  "Study  of  the  Interaction  of  Light 
with  Rough  Metal  Surfaces,"  Phys.  Rev.  B  2^:309-329  (1970). 

Beard,  C.  I.,  "Coherent  and  Incoherent  Scattering  of  Microwaves 
from  the  Ocean,"  IRE  Trans.  AP9:470-483  (1961). 

Beckmann,  P.,  "A  New  Approach  to  the  Problem  of  Reflection  from  a 
Rough  Surface,"  Acta.  Tech.  CSAV.  2^: 311-355  (1957). 

Bennett,  H.  E.,  "Specular  Reflectance  of  Aluminized  Ground  Glass  and 
the  Height  Distribution  of  Surface  Irregularities,"  J.  Opt.  Soc. 
Am.  53:1389-94  (1963). 

Bennett,  H.  E.,  and  J.  M.  Bennett,  "Precision  Measurements  in  Thin 
Film  Optics,"  in  Physics  of  Thin  Films,  G.  Hass  and  R.  E.  Thun, 
eds.,  Academic  Press,  N.  Y.,  Vol.  £,  pp.  1-96  (1967). 

Bennett,  H.  E.,  and  J.  0.  Porteus,  "Relation  between  Surface  Roughness 
and  Specular  Reflectance  at  Normal  Incidtnce,"  J.  Opt.  Soc.  Am. 

51_:  123-9  (1961). 

Bennett,  J.,  "Use  of  Interferometry  for  Determining  the  rms  Roughness, 
Autocovariance  Function,  and  Other  Statistical  Properties  of 
Optical  Surfaces,"  (Abstract)  J.  Opt.  Soc.  Am.  64:1362  (1974). 


155 


Blake,  L.  V.,  "Reflection  of  Radio  Waves  from  a  Rough  Sea,"  Proc. 

IRE  38:301-304  (1950) . 

Bloembergen,  N.,  and  Y.  R.  Shen,  "Coupling  of  Light  with  Phonons, 
Magnons,  and  Plasmons,"  Phys.  of  Quantum  Elect.  Conference  Proc., 

P.  L.  Kelley,  et  al.,  eds.,  San  Juan,  Puerto  Rico,  McGraw-Hill, 

N.  Y.,  pp.  119-128  (1965). 

Bojarski,  N.  N.,  "Electromagnetic  Inverse  Scattering,"  Mathematics 
of  Profile  Inversion  Workshop,  Ames  Research  Center,  Moffett 
Field,  California,  L.  Colin,  ed.,  pp.  6-34  to  6-37,  NASA  TM  X-62, 
150,  Available  from  the  National  Technical  Information  Service, 

No.  N73-11615,  (1971). 

Bracewell,  R.  N.,  The  Fourier  Transform  and  its  Applications, 
McGraw-Hill,  N.  Y.  (1965). 

Breault,  R.  P.,  and  B.  B.  Fannin,  Optical  Sciences  Center,  University 
of  Arizona,  Tucson,  Private  Communication  (1973) . 

Chenmoganadara,  T.  K.,  "On  the  Specular  Reflection  from  Rough  Surfaces," 
Phys.  Rev.  13:96-101  (1919). 

Church,  E.  L.,  and  J.  M.  Zavada,  "Residual  Surface  Roughness  of 
Diamond -Turned  Optics,"  Appl.  Opt.  1£:1788  (1975). 

Crowell,  J.,  and  R.  H.  Ritchie,  "Surface  Plasmon  Effect  in  the 
Reflectance  of  a  Metal,"  J.  Opt.  Soc.  Am.  60:974-799  (1970). 

Curcio,  M.  E.,  "Evaluation  of  Low  Scatter  Technology  for  Aspheric 
Metal  Mirrors,"  Oral  Presentation  at  SPIE's  19th  Annual  Inter¬ 
national  Technical  Symposium,  August  1975,  San  Diego,  CA  (1975). 

Daniels,  F.  B.,  "A  Theory  of  Radar  Reflection  from  the  Moon  and 
Planets,"  J.  Geophys.  Res.  66:1781-1788  (1961). 

Daud6,  A.,  A.  Savary,  and  S.  Robin,  "Effects  of  Different  Rough¬ 
nesses  on  the  Excitation  by  Photons  of  the  Surface  Plasmon 
of  Aluminum,"  J.  Opt.  Soc.  Am.  62:1-5  (1972). 

Davies,  H. ,  "The  Reflection  of  Electromagnetic  Waves  from  a  Rough 
Surface,"  Proc.  IEEE  10j^:209-14  (1954). 

DeBell,  M.  A.,  and  J.  E.  Harvey,  "Recent  Developments  in  Surface 
Scatter  Studies,"  Abstract,  J.  Opt.  Soc.  Am  64:1404  (1974). 

Decker,  D.  L.,  H.  E.  Bennett,  and  J.  M.  Bennett,  "Light  Scattering 
and  Absorption  in  Infrared  Laser  Mirrors,"  Oral  Presentation, 

SPIE's  19th  Annual  International  Technical  Symposium,  August, 

1975,  San  Diego,  California  (1975). 


156 


Eckart,  C.,  "The  Scattering  of  Sound  from  the  Sea  Surface," 

J.  Acoust.  Soc.  Am.  £5:566-570  (1953). 

Elson,  J.  M.,  and  R.  H.  Ritchie,  "Photon  Interactions  at  a  Rough 
Metal  Surface,"  Phys.  Rev.  B  £:4129-4138  (1971). 

Evans,  J.  V.,  and  G.  H.  Pettengill,  "The  Scattering  Behavior  of 
the  Moon  at  Wavelengths  of  3.6,  68,  and  784  Centimeters," 

J.  Geophys.  Res.  68^: 423-447  (1963). 

Feinstein,  J.,  "Some  Stochastic  Problems  in  Wave  Propagation," 

Part  I.  Trans.  IRE  AP  2: 23-30  (1954). 

Fuks,  I.  M.,  "Theory  of  Radio  Wave  Scattering  at  a  Rough  Sea  Surface," 
Soviet  Radiophysics  9_:  513-519  (1966). 

Fung,  A.  K.,  "Character  of  Wave  Depolarization  by  a  Perfectly  Con¬ 
ducting  Rough  Surface  and  Its  Application  to  Earth  and  Moon 
Experiments,"  Planetary  and  Space  Science  £5:1337-3147  (1967). 

Fung,  A.  K.,  and  R.  K.  Moore,  "Effects  of  Structure  Size  on  Moon 
and  Earth  Radar  Returns  at  Various  Angles,"  J.  Geophys.  Res. 

69: 1075-1081  (1964). 

Goodman,  J.  W.,  Introduction  to  Fourier  Optics,  McGraw-Hill,  San 
Francisco,  CA  (1968). 

Hagfors,  T.,  "Backscattering  from  an  Undulating  Surface  with  Appli¬ 
cations  to  Radar  Returns  from  the  Moon,"  J.  Geophys,.  Res. 
69:3779-3784  (1964). 

Hopkins,  H.  H.,  Wave  Theory  of  Aberrations,  Oxford  Clarendon  Press, 
London  (1950) . 

Isakovich,  M.  A.,  "The  Scattering  of  Waves  from  a  Statistically 

Rough  Surface,"  Zhurn.  Eksp.  Theor.  Fiz.  23:305-314,  Translation 
available  from  Defense  Scientific  Information  Service,  DRB 
Canada,  T  112  R  (1952). 

Judd,  D.  B.,  "Terms,  Definitions,  and  Symbols  in  Reflectometry," 

J.  Opt.  Soc.  Am.  57:445-452  (1967).' 

Lalor,  E.,  "Conditions  for  the  Validity  of  the  Angular  Spectrum  of 
Plane  Waves,"  J.  Opt.  Soc.  Am.  £8:1235  (1968). 

Leinert,  C.,  and  D.  Kliipelberg,  "Stray  Light  Suppression  in  Optical 
Space  Experiments,"  Appl.  Opt.  £3:556  (1974). 


McKenney,  D.  B.,  G.  R.  Orme,  and  L.  P.  Mott,  "Light  Scattering  by 
Thin  Film  Coatings,"  Final  Report,  Optical  Sciences  Center, 
University  of  Arizona,  Tucson  (1972) . 

Moffitt,  F.  H.,  Photogrammetry,  International  Textbook  Co.  (1959). 

Mott,  L.  P.,  The  Effect  of  Surface  Roughness  on  the  Optical  Proper¬ 
ties  of  All-dielectric  Interference  Filters,  Master's  thesis. 
University  of  Arizona,  Tucson  (1971). 

Nankivell,  J.  F.,  "The  Theory  of  Electron  Stereo  Microscopy,"  Optic 
20j_  171  (1963) . 

Nicodemus,  F.  E.,  "Reflectance  Nomenclature  and  Directional  Reflec¬ 
tance  and  Emissivity,"  Appl.  Opt.  j):  1474-5  (1970). 

Orme,  G.  R.,  "Measurement  of  Small  Angle  Scatcer  from  Smooth  Sur¬ 
faces,"  Optical  Sciences  Center  Tech.  Rpt.  7£,  University  of 
Arizona,  Tucson  (1972) . 

Papoulis,  A.,  Probability,  Random  Variables,  and  Stochastic  Processes, 
McGraw-Hill,  N.  Y.  (1965). 

Peake,  W.  H.,  "Interaction  of  Electromagnetic  Waves  with  Some  Natural 
Surfaces,"  IRE  Trans.  AP-7  (1959). 

Perkin-Elmer  Final  Report  pref  red  for  George  C.  Marshall  Space 
Flight  Center  under  Contract  NAS8-30638,  "Study  of  Ultraviolet 
and  Visible  Scattered  Light  Effects  on  the  Optical  Performance 
of  the  Large  Space  Telescope,"  (1975). 

Porteus,  J.  0.,  "Relation  between  the  Height  Distribution  of  a  Rough 
Surface  and  the  Reflectance  at  Normal  Incidence,"  J.  Opt.  Soc. 

Am.  53:1394-1402  (J963). 

Ratcliffe,  J.  A.,  "Some  Aspects  of  Diffraction  Theory  and  Their 
Application  to  the  Ionosphere,"  in  A.  C.  Strickland,  ed., 

Reports  on  Progress  in  Physics,  Vol.  XIX,  The  Physical  Society, 
London  (1956). 

Rayleigh,  Lord,  "Polish,"  Nature  6£: 385-388  (1901). 

Rayleigh,  Lord,  The  Theory  of  Sound,  Vol.  II,  Dover,  N.  Y.  (1945). 

Scheele,  S.  R.,  Scattering  Characteristics  of  Mirrors  and  the 
Associated  Inverse  Scattering  Problem,  Ph.  D.  dissertation. 
University  of  California  at  Los  Angeles  (1973). 


158 


Shack,  R.  V.,  "Interaction  of  an  Optical  System  w:th  the  Incoming 
Wavefront  in  the  Presence  of  Atmospheric  Turbulence,"  Optical 
Sciences  Center  Tech.  Rpt.  1£,  University  of  Arizona,  Tucson 
(1967). 

Shack,  R.  V.,  and  M.  A.  DeBell,  "Surface  Scatter  Study,"  Optical 
Sciences  Center  final  report  prepared  for  the  Space  and  Missile 
Systems  Organization  under  Contract  F04701-72-C-0181 ,  University 
of  Arizona,  Tucson  (1974) . 

Shack,  R=  V.,  and  J.  E.  Harvey,  "An  Investigation  of  the  Distribution 
of  Radiation  Scattered  by  Optical  Surfaces,"  Optical  Sciences 
Center  final  report  prepared  for  the  Space  and  Missile  Systems 
Organization  under  Contract  F04701-74-C-0083,  University  of 
Arizona,  Tucson  (1975). 

Shevlin,  C.  M.,  Department  of  Metallurgy,  University  of  Arizona, 
Tucson,  Private  Communication  (1974). 

Spetner,  L.  M.,  "A  Statistical  Model  for  Forward  Scattering  of  Waves 
of  a  Rough  Surface,"  Trans.  IRE  AP  6^: 88- 94  (1958). 

Stover,  J.  C.,  "Roughness  Characterization  of  Smooth  Machined  Surfaces 
by  Light  Scattering,"  Appl.  Opt.  1_4:1796  (1975). 

Thompson,  S.  P.,  Treatise  on  Light,  Macmillan  and  Co.,  London  (1912). 

Twersky,  V.,  "On  Scattering  and  Reflection  of  Electromagnetic  Waves 
by  Rough  Surfaces,"  TRE  Trans.,  pp.  81-90  (1957). 

Vachaspati,  "Scattering  of  Light  by  Free  Electrons,  Intensity 

Dependent  Phase  Shift  and  Frequency  Change,"  Indian  J.  Pure  and 
Appl.  Phys.  2:373-376  (1964). 

Weyl,  H.,  "Ausbreitung  elektromagnetischer  wellen  uber  einen 
ebenen  leiter,"  Ann.  Phys.  60:481-500  (1919). 

Young,  R.  P.,  "Metal  Optics  Scatter  Measurements,"  Oral  Presentation 
at  SPIE's  19th  Annual  International  Technical  Symposium,  San 
Diego,  CA,,  August  (1975). 


